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Abstract. By means of three-dimensional hydrodynamic 
simulations with a Eulerian PPM code we investigate the 
formation and the properties of the accretion torus around 
the stellar mass black hole which originates from the merg- 
ing of two neutron stars. The simulations are performed 
with four nested cartesian grids which allow for both a 
good resolution near the central black hole and a large 
computational volume. They include the use of a physical 
equation of state as well as the neutrino emission from 
the hot matter of the torus. The gravity of the black hole 
is described with a Newtonian and alternatively with a 
Paczyhski-Wiita potential. In a post-processing step, we 
evaluate our models for the energy deposition by vv an- 
nihilation around the accretion torus. 

We find that the torus has a mass between several 
1(T 2 M© and a few 1O _1 M with maximum densities 
around 10 12 g cm -3 and maximum temperatures of about 
lOMeV (entropies around 5 fee per nucleon). Correspond- 
ingly, the neutrino emission is huge with a total luminos- 
ity near 10 53 ergs _1 . Neutrino-antineutrino annihilation 
deposits energy in the vicinity of the torus at a rate of (3- 
5) x 10 50 ergs . It is most efficient near the rotation axis 
where 10 to 30% of this energy or up to a total of 10 49 erg 
are dumped within an estimated emission period of 0.02- 
0.1 s in a region with a low integral baryonic mass of about 
10~ 5 Mq. This baryon pollution is still dangerously high, 
and the estimated maximum relativistic Lorentz factors 
r — 1 are around unity. The conversion of neutrino en- 
ergy into a pair plasma, however, is sufficiently powerful to 
blow out the baryons along the axis so that a clean funnel 
should be produced within only milliseconds. Our models 
show that vv annihilation can yield the energy to account 
for weak, short gamma-ray bursts, if moderate beaming 
is involved. In fact, the barrier of the dense baryonic gas 
of the torus suggests that the low-density e ± 7 plasma is 
beamed as axial jets into a fraction fa = 2<5f2/(4-7r) be- 
tween 1/100 and 1/10 of the sky, corresponding to open- 
ing half-angles of roughly ten to several tens of degrees. 
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Thus 7-burst energies of E 1 w E uP /fa <> 10 5U -10 51 erg 
seem within the reach of our models (if the source is 
interpreted as radiating isotropically), corresponding to 
luminosities around 10 51 ergs _1 for typical burst dura- 
tions of 0.1-1 s. Gravitational capture of radiation by the 
black hole, redshift and ray bending do not reduce the jet 
energy significantly, because most of the neutrino emis- 
sion comes from parts of the torus at distances of several 
Schwarzschild radii from the black hole. Effects associated 
with the Kerr character of the rapidly rotating black hole, 
however, could increase the 7-burst energy considerably, 
and effects due to magnetic fields might even be required 
to get the energies for long complex gamma-ray bursts. 

Key words: gamma rays: bursts - elementary parti- 
cles: neutrinos - stars: neutron - binaries: close - hydro- 
dynamics 



1. Introduction 



In a sequence of prec eding papers (Ruffert et al 



1996 



1997| ; Ruffert & Janka |l998a| , |l998b|) we have shown that 



the neutrino emission associated with the dynamical phase 
of the merging or collision of two neutron stars is power- 
ful, but too short to provide the energy for gamma-ray 
bursts by neutrino-antineutrino annihilation. Significant 
heating of the coalescing stars occurs only after they have 
plunged into each other, and the neutrino luminosities can 
rise to several 10 53 ergs _1 and even exceed 10 54 ergs _1 
in case of the more violent collisions. After a few mil- 
liseconds, however, the compact massive remnant of the 
merger will most likely collapse to a black hole. If that did 
not happen, the remnant's continuing neutrino emission 
would drive a dense baryonic wind off its surface which 
would lead to a sizable mass loss but non-relativistic ex- 
pansion (Woosley & Baron |1992| Woosley p"993b| , Hernanz 
et al. 1994, Qian & Woosley 1996), a situation which is not 
favorable for producing gamma-ray bursts which require 
relativistic Lorentz factors T *b 100 (Paczyhski 199C). 

If the merger remnant collapses to a black hole, some 
matter remains in an accretion disk or torus around the 
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black hole. Our hydrodynamic simulations (Ruffert et al. 
1996 , Ruffert & Janka 1998b ) have given hints that about 



0.1 Mq of matter might obtain enough angular momentum 
during the merging of the neutron stars to resist immedi- 
ate collapse into the black hole. In this case a funnel with 
low baryon density can develop along the system axis. 
On the other hand, the large angular momentum ensures 
that the torus matter is swallowed by the black hole on 
a time scale much longer than the dynamical time scale. 
Therefore there could be enough time for this material 
to radiate away a fair fraction of its gravitational bind- 
ing energy in neutrinos, even if the densities become so 
high that neutrinos get trapped and can escape only on a 
diffusion time scale. A similar situation could result from 
the merging of a neutron star with a black hole (Lee & 
Kluzniak |1995j |1998| ; Kluzniak & Lee |1998j ; Eberl |1998a ; 
Eberl et al. 1998b| ), from the collaps e of a very m assive , 



rapidly rotating star ( Woosle y |1993a[ Popham et al. |1998| , 
MacFadyen & Woosley |l998|) , or from the coalescence of 
a neutro n star/black hole with a white dwarf (Fryer et 
al. 1998 ) or with the helium core of its red giant compan- 
ion (Fryer & Woosley 1998] ). 



All these events are estimated to occur at rates which 
can account for the observed frequency of gamma-ray 
bursts (about one burst per day). Also, the huge amount 
of gravitational binding energy released during the accre- 
tion process of up to several solar masses of gas into the 
black hole is hoped to be able to explain the energetics 
of even the most distant cosmological gamma-ray bursts 
(e.g., GRB981214, see Kulkarni et al. |1998| ). Moreover, the 
compactness of the stellar-mass black hole could naturally 
produce the rapid variability on time scales of milliseconds 
observed in many bursts. For these reasons, massive ac- 
cretion disks or thick accretion tori around stellar-mass 
black holes are considered as possible cosmological origin 
of the enigmatic gamma-ray bursts, powered by neutrino- 
antineutrino annihilation or by magnetically driven energy 
rel ease ( e.g., Paczyhski 1986 ; Goodman 1986 ; Goodman et 
al.|l987| ; Eichler et al. |1989|; P aczyhski [L991[ N arayan et al. 
1992 ; Mesz aros & Rees 1993| ; Wo osley 19934 Jaroszyhski 



1992, 1996; Mochkovitch et al. 1993, 1993; Thompson 



1994; Witt et al. 1994; Janka fc Ruffert 1996; Meszaros 



& Rees |1997t Popham et al. |1998j ; Meszaros et al. |1998|). 

In this paper we simulate the formation of the accre- 
tion torus after two neutron stars have merged, and as- 
sume that the compact remnant with a baryonic mass of 
about 3 Mq has collapsed into a black hole. The initial 
model is taken from our merger simulations (Ruffert & 
Janka 1998b ) where the central, massive object is replaced 
by a vacuum boundary at a radius equal to twice the 
Schwarzschild radius of the mass dumped into the black 
hole. We follow the evolution of the left-over material in 
the surroundings of the black hole until it either has set- 
tled into the disk or has been swallowed by the black hole. 
At the end of our three-dimensional computations, the 
torus has reached a quasi-stationary state and its further 



evolution is governed by the viscous transport of angu- 
lar momentum which depends on the uncertain value of 
the disk viscosity. Kerr effects associated with the rotation 
of the black hole and magnetic fields are not taken into 
account in our models. An extensive, general investiga- 
tion of relativistic steady-state accretion from tori around 
hyper-accreting stellar black holes for different accretion 
rates and disk viscosities in Schwarzschild and Kerr geom- 
etry was recently published by Popham et al. (1998). Our 
simulations focus on the situation that emerges from the 
mergings of compact binary systems of neutron stars and 
black holes, and they are intended to help constrain the 
large parameter space and to yield insight into the torus 
properties and non-stationary aspects of the evolution. 

In particular, our simulations aim at answering the 
following questions: How much mass remains in the ac- 
cretion torus? What is the relativistic rotation parameter 
a = Jc/(GM 2 ) of the black hole? What are the properties 
of the accretion torus, its density, temperature, neutrino 
luminosity? How much mass pollutes the surroundings 
of the accretion torus, in particular, does an effectively 
baryon-free funnel form along the system axis? How ef- 
ficient is neutrino-antineutrino annihilation in depositing 
energy in the regions with low baryon density? Can we 
make estimates of the mass accretion rate into the black 
hole and the corresponding lifetime of the torus? What 
are the implications for producing gamma-ray bursts by 
neutrino-antineutrino annihilation? Is pair-plasma ejected 
in a jet and how large will its opening angle be? Is there 
enough variability of the energy release at the central 
source to account for the observed time structure of the 
gamma-ray burst light curves? 

The paper is organized as follows. In Sect. ||the compu- 
tational procedures are summarized which are used in our 
simulations of the torus formation after neutron star merg- 
ing. The initial model for these simulations is briefly de- 
scribed and the different investigated cases are introduced 
with their distinguishing parameters. Section || contains a 
description of the dynamical evolution of the torus from 
the beginning of the simulations until a quasi-stationary 
state was reached. In Sect. || the properties of the accretion 
tori at the end of the computations are described. The re- 
sults on the neutrino emission are presented in Sect. p| and 
those for neutrino-antineutrino annihilation in Sect. |6J. In 
Sect. ^ the hydrodynamic results are used to estimate the 
numerical viscosity which determines our torus models; 
the general relativistic effects in the neutrino-antineutrino 
annihilation are discussed as well as the importance of 
neutrino-electron/positron scattering for the heating of 
the pair-plasma cloud that is formed by vv annihilation. 
Section || concludes the paper with a summary and a dis- 
cussion of the implications of our results for gamma-ray 
burst scenarios involving massive accretion tori around 
stellar mass black holes. 
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Table 1. Some parameters for the torus evolution models with Newtonian potential (Models Bl, B2, B4, and BIO) and 
with Paczyhski-Wiita potential (Models Bl, B2, and BIO), resp ectivel y. The dynamical simulations were started with 
initial conditions as given from Model B64 of Ruffert & Janka ( 1998b|) at different times i; n ;t and the evolutions were 



followed over time intervals Ai ca i. All other quantities refer to the conditions found at the end of the simulations. M p< n 
is the mass on the grid with density less than 10 11 g cm -3 , M v the total gas mass on the computational grid, and Md the 
mass of the gas with specific angular momentum larger than the Kepler limit jjj = v£J epler (3i? s )3i? s = y/EGM/c in case 
of the Newtonian potential or jp W = v^^ )lcl (3R s )3R s = §jj!j for the Paczyhski-Wiita potential; the Kepler velocities 
were evaluated from Eqs. (||) and @ using for M the sum of black hole and total gas mass on the grid at the end of 
the simulation. !T max is the maximum gas temperature in energy units. L Ve denotes the electron neutrino luminosity 
near the end of the simulation, L Pc the corresponding electron antineutrino luminosity, and L Vi the luminosity of each 
individual flavor of heavy- lepton neutrino (z/^, D^, v r or P T ). The sum of all individual neutrino luminosities is given 
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2. Computational procedures, initial conditions 
and different models 



In Ruffert & Janka ( |l998bj ) we have calculated a series 
of neutron stars merger models with varied neutron star 
masses, neutron star mass ratios, neutron star spins, and 
initial conditions (temperature, entropy) in the coalescing 
stars. In all of these simulations the post-merging con- 
figuration consisted of a compact central object with a 
mass of about 3 Mq and a typical density of the order of 
10 14 gcm -3 which was surrounded by an extended cloud 
of more dilute gas, having a mass of a few 0.1 Mq and a 
characteristic mean density around 10 11 -10 12 g cm~ 3 . To 
estimate the gas mass which might be able to stay in an 
accretion torus after the massive object has collapsed into 
a black hole, we compared the specific angular momentum 
j of the matter to the Keplerian angular momentum j* for 
a test particle with non-zero mass that orbits around the 
Schwarzschild black hole on the last stable circular orbit 
at 3 Schwarzschild radii, 3i? s = 6GM/c 2 . Taking for sim- 
plicity M as the total (gas) mass on the grid, we found 
that between several 10~ 2 Mq and a few 10 _1 Mq ful- 
fill this criterion : j > j* = v£ cplcr (3i? s )3i? s = V6GM/c 
where V^epler i s the Keplerian velocity and the superscript 
(or subscript) N indicates the use of a Newtonian gravi- 
tational potential. The amount of mass which has suffi- 
ciently large angular momentum to resist immediate (i.e., 
on a dynamical time scale) accretion into the black hole 
depends on the neutron star masses and mass ratio as well 
as on the neutron star spins which contribute to the total 
angular momentum of the binary system. For the simu- 



lations of the formation of the accretion torus which we 
describe in this paper, the neutron star merger Model B64 
of Ruffert & Janka (1998b ) was used as an initial condi- 
tion (defined by the distributions of density, temperature 
and electron fraction, and by the velocity field at a certain 
chosen time during the evolution of Model B64). Due to 
the assumed corotation of the two 1.6 Mq (baryonic mass) 
neutron stars before merging, the angular momentum was 
largest in this model and correspondingly, the estimated 
possible torus mass was maximal. 

The three-dimensional computations of neutron star 
mergings were performed with a Newtonian hydrodynam- 
ics code based on the Piecewise Parabolic Method (PPM) 
of Colella & Woodward (1984) with at least four levels of 
nested grids (Ruffert 1992) to ensure both high resolution 
at the neutron stars and a large computational volume. 
The code includes the effects of gravitational-wave emis- 
sion and their back-reaction on the hydrodynamic flow ac- 
cording to Blanchet et al. ( 199C| ) (see Ruffert et al. 1996). 
In addition, we implemented a calibrated neutrino leakage 
scheme (Ruffert et al. 1996) in order to calculate the en- 
ergy and lepton number loss by neutrino emission from the 
heated neutron star matter. The latter is described by the 
finite-temp eratu re nuclear equation of state of Lattimcr 
& Swesty C tL99l|) using the Skl80 nuclear force parameter 
set (Swesty et al. |l994[) . 

The torus simulations presented in this paper are done 
with the same code and the same input physics. Because 
of the spherical symmetry of the mass in the black hole 
and the relatively small mass of the accretion torus, time 
derivatives of the quadrupole moment are small and grav- 
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Fig. 1. Maximum absolute value of the gravitational po- 
tential produced by the mass distribution of Model B64 of 



Ruffert & Janka (1998b) as function of time, normalized 
to c 2 . This parameter, multiplied by 2, is a measure of the 
importance of relativistic corrections, and is very small in 
the limit of Newtonian gravity. The time is measured from 
the beginning of the simulation of the neutron star merger 
Model B64. 

itational-wave production does not play an important role. 
Although taken into account in our simulations, we shall 
therefore not report data of the gravitational-wave emis- 
sion here. In one set of our models (Bl, B2, B4, and 
BIO) the black hole potential, which dominates the grav- 
itational field at the torus (whose self-gravity is only a 
minor contribution), is represented by a Newtonian po- 
tential, 



GM] 



BH 



(1) 



In a second sequence of models (Bl, B2, and BIO) the grav- 
itational potential is described by the Paczyhski-Wiita ex- 
pression, 



GM BH 
r-IL 



(2) 



(Paczyhski & Wiita |1980|) . This allows one to reproduce 
the existence and the effects of a last stable circular orbit 
at a radius of 3i? s = 6GA/bh/c 2 where the specific angu- 
lar momentum j PW = ™Ke^ler( r ) = V GA/bh^ 3 / (r - R s ) 2 
has an absolute minimum. We hope that this approxima- 
tion, although crude, can give us some indication of the 
sensitivity of our results to the inclusion of proper general 
relativity in the modeling. 

The different models of each set, B# and 23#, respec- 
tively, are discerned by the different times t init (in mil- 
liseconds given by the numbers in the model names) at 
which the compact object that formed during the merg- 
ing of the two neutron stars in Model B64 of Ruffert & 



Janka ( 1998b ) was removed and replaced by a gravitat- 
ing "vacuum sphere" . The times ii n it are measured from 
the start of the simulation of Model B64 and are listed 
in Table [fl. If a black hole forms from the compact cen- 
tral body of the merger remnant with a baryonic mass 
of nearly 3 M Q , we expect this to happen at about the 
time when the most compact state is reached and thus 
the gravitational potential becomes strongest. If support 
by centrifugal forces or thermal pressure played a role, the 
collapse might be delayed. Figure [l] shows that the value 
of the parameter $n/c 2 — twice of which gives a rough 
measure for the importance of general relativistic gravity 
— plateaus at about 2 ms after the simulation of the merg- 
ing of the neutron stars was started. Therefore we chose 
the earliest moment of the possible black hole formation 
and the onset of our torus computations at around this 
time, and investigated also models for later collapse times 
at about 4 ms and 10 ms, respectively. The model runs 
were continued until the accretion rate into the black hole 
had reached such a low value that further changes of the 
torus properties would happen over a much longer period 
than the dynamical time scale of the system. The com- 
puted evolution times Ai ca i of all models are also given in 
Table |l|. The subsequent quasi-stationary evolution pro- 
ceeds on the time scale of viscous transport of angular 
momentum which depends on the unknown value of the 
disk viscosity and in general is too long to be followed by 
our three-dimensional, hydrodynamic simulations with an 
explicit code. 

The inner vacuum sphere that represents the black hole 
at the center of the computational grid is set to a radius 
of 2R S . The Schwarzschild radius R s is initially computed 
for 90% of the total gas mass on the grid and the mass of 
the gas inside this radius is collected into the black hole to 
determine its gravitational potential. During the following 
evolution, the black hole mass, momentum and angular 
momentum are updated by adding the corresponding val- 
ues of the matter which is advected through the sphere 
at 2R S . From the current value of the black hole mass, 
the new radius of the vacuum sphere and the new grav- 
itational potential of the black hole are calculated. The 
loss of mass, momentum, angular momentum and energy 
from the gas outside the black hole boundary are also 
monitored during the simulations. In the grid zones that 
are located inside the vacuum sphere (these zones are not 
removed from the hydrodynamic grid), the mass density is 
continuously reset to a negligibly small but finite value of 
10 8 g cm -3 and a correspondingly very small value of the 
pressure is used. The decision to put the vacuum bound- 
ary at the radius 2R S was influenced by the facts that on 
the one hand the gravitational potential in the Paczyhski- 
Wiita case diverges when r decreases towards 1R S , and 
that on the other hand the gas velocities come close to the 
speed of light already near 2R S and therefore the nonrela- 
tivistic treatment of the hydrodynamics can definitely not 
be applied any longer. Moreover, the Courant-Friedrich- 
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Fig. 2. Distribution of mass and angular momentum in 
the remnant of the neutron star merger of Model B64 of 
Ruffert & Janka (1998b) at time t = 10ms. This defines 
the initial condition used for some of the torus simulations 
presented here. The upper thin solid line gives the cumu- 
lative gas mass (in solar masses), the lower thin solid line 
the cumulative angular momentum (in 10 49 g cm 2 s^ 1 ), 
and the thick solid line the relativistic rotation param- 
eter a(r) = J(r)c/(GMg as (r)) as functions of radius r. 
(M gas (r) is the gas mass inside r and J(r) the correspond- 
ing angular momentum perpendicular to the equatorial 
plane.) 
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Fig. 3. Mass accretion rates of the black hole (in 
MqHis -1 ) as functions of time for all models listed in 
Table |l|. The thick lines which start at times t = 1.84 ms, 
2.60 ms, 4.09 ms, and 10.0 ms correspond to the Newto- 
nian Models Bl, B2, B4, and B10, respectively, and the 
thin lines starting at times t = 1.84ms, 2.60ms, and 
10.0 ms correspond to Models Bl, B2, and 610, respec- 
tively. 
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Fig. 4. The mass of the central black hole as function of 
time for Model B10 (bold line) and Model B10 (thin line), 
showing an increase due to accretion of surrounding gas. 



Fig. 5. Typical time scales for the changes of black hole 
mass (solid lines) and torus mass (dashed lines) as func- 
tions of time for Model B10 (bold lines) and Model B10 
(thin lines). 
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Fig. 6. Different masses as functions of time for the Newtonian Models Bl, B2, B4, and B10 (left plot) and for 
the different Paczyhski-Wiita Models Bl, B2, and B10 (right plot) (distinguishable by their start and stop times, see 
Table [|). M v (bold solid lines) denotes the gas mass on the grid, M u (dashed lines) is the cumulative mass which 
leaves the grid during the simulation and becomes unbound, i.e., which fulfils the criterion that its total specific energy 
as the sum of its specific internal, kinetic, and potential energies is positive. The dotted curves with labels Mio, Mu, 
and M12 show the gas masses with densities below 10 10 gcm -3 , 10 1:L gcm _3 , and 10 12 gcm -3 , respectively. In case of 
the Newtonian potential (left plot), the thin solid lines labeled with represent the mass on the grid with specific 
angular momentum larger than the Kepler value at 3 Schwarzschild radii, i.e., j > = \^6GM/c, where for M the 
sum of the black hole mass and gas mass on the grid was used at all times. 
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Fig. 7. Maximum density on the grid as function of time 
for the Newtonian Models B2, B4, and B10 (bold lines) 
and for the Paczyhski-Wiita Models B2 and B10 (thin 
lines). 



Fig. 8. Maximum temperature on the grid as function of 
time for the Newtonian Models B2, B4, and B10 (bold 
lines) and for the Paczyhski-Wiita Models B2 and B10 
(thin lines). 



Lewy timestep of the explicit computation is limited by 
the small values in this region on the finest grid. From a 
physics point of view, the choice of the black hole bound- 
ary at 2i? s can be justified because the generated black 
hole is not of extreme Kerr type but its relativistic rota- 
tion parameter a is of order 0.4-0.5 (Fig. |^). In this case 
the inner edge of the torus should be located at just around 



two times the horizon radius (see Usui et al. 1998). Interior 
to this radius the orbits are unstable and the infall veloc- 
ities become very large, i.e., the matter moves essentially 
radially inward. Our Newtonian code does not take into 
account effects due to the rotation of the black hole on the 
surrounding space and matter, e.g., frame dragging or the 
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Fig. 9. Average value of the specific internal energy of Fig. 10. Total kinetic, thermal, and gravitational poten- 
the matter on the grid as function of time for Model BIO tial energies as functions of time for Model BIO (bold 



(bold line) and Model 610 (thin line). 



line) and Model BIO (thin line). 



dependence of the radius of the innermost stable orbit on 
the rotation parameter of the black hole. 

The simulations are performed with four levels of 
nested cartesian grids, each having 64 zones per dimen- 
sion in the orbital plane and 16 perpendicular to the or- 
bital plane (here we make use of equatorial symmetry 
and of the smaller extension of the torus in the verti- 
cal direction). A zone on the finest grid has a length of 
Ax = Ay = Az = 0.64 km, and the size of the largest grid 
is 328 km x 328 km x 82 km. 



3. Dynamical evolution 

Figures ^[-jlfj| show the evolution of the remnant of the 
merged neutron stars after the compact central object has 
been replaced by a vacuum sphere to represent the black 
hole. The different simulations, discriminated by the use 
of a Newtonian or Paczyhski-Wiita type of gravitational 
potential and by different times of the assumed black hole 
formation, are represented by different line styles. Time is 
measured from the start of the simulation of Model B64 



of Ruffert & Janka fll998b| ). 

In Fig. H the mass accretion rates of the black hole 
are given as functions of time for all models, in Fig. |] the 
evolution of the black hole mass is displayed for the New- 
tonian and Paczyhski-Wiita models which were started at 
t = 10ms (Models BIO and BIO, respectively), and in 
Fig. U the corresponding time scales for the changes of 
black hole and torus masses are shown. Initially, the black 
hole swallows the surrounding mass at rates as high as 
several solar masses per millisecond, but within a dynam- 
ical time scale of only about 1 ms the accretion rates set- 



tle to much lower values between 5 M© s _1 and 7M© s _1 . 
The peak rates as well as the rates towards the end of 
the simulations are similar in all models. The amount of 
surrounding gas which is dynamically accreted into the 
black hole is larger for the Paczyhski-Wiita case, and the 
black hole mass grows correspondingly faster (Fig. ^). 
When the simulations are stopped after a quasi-stationary 
situation has been reached, the torus is therefore about 
10 times more massive in the Newtonian models (torus 
mass w 0.25-0.35 M©; compare Table 0) than in 

those with Paczyhski-Wiita potential (M^ w w 0.005- 
0.035 M©; compare Table [|). Note, however, that the sim- 
ulations with the Paczyhski-Wiita black hole potential 
might well underestimate the torus mass. Since these sim- 
ulations were started with an initial model that resulted 
from a Newtonian computation of neutron star merging, 
the initial angular momentum of the gas may be lower 
than would have been obtained in a simulation with the 
stronger relativistic potential. Towards the end of the sim- 
ulations, the time scales of the changes of black hole mass 
and accretion torus mass level off at values much larger 
than the dynamical time scale. The accretion time scale 
of the Newtonian torus grows to about 60 ms, whereas 
the accretion time scale for the less massive torus in the 
Paczyhski-Wiita case is approximately one order of mag- 
nitude shorter (Fig. ||). 

Figure || demonstrates that the torus mass — given 
as the gas mass M v on the grid — approaches nearly the 
same value in the quasi-stationary state when the moment 
of black hole formation is ii n j t > 2 ms. In contrast, if the 
black hole is assumed to form early after the merging of the 
neutron stars (ti n it < 2ms; Models Bl and Bl), it swal- 
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Fig. 11. Structure of the Newtonian Model B10 in the equatorial plane near the end of the simulation. Panel a 
shows contours of the baryonic mass density (measured in gem -3 ) together with the flow field indicated by velocity 
vectors, panel b displays contours of the temperature (in MeV), panel c of the electron fraction Y e , and panel d of 
the entropy per nucleon. The length scale of the ve locity vectors and the time elapsed since the beginning of the 
simulation of Model B64 of Ruffert & Janka ( 1998b ) are given in the insert in the top right corner of the figures. 
The density contours are spaced logarithmically with intervals of 0.5 dex (bold: log(p [gem -3 ]) = 10.0, 11.0, 12.0), 
the temperature contours are spaced linearly, starting with 1 MeV and 2 MeV, and then continuing in steps of 2 MeV 
(bold: T = 10 MeV), the Y e contours are spaced linearly with intervals of 0.02 (bold: 0.02, 0.06, 0.10, 0.16), and the 
entropy levels were chosen to be 1, 2, 3, 4, 5, 6, 8, 10, 12, 14, 16, and 20 (bold: 6, 10, 20). Some of the contours are 
labeled with their respective values, and darker grey shading indicates higher values of a quantity. The circle around 
the center marks the inner "vacuum" boundary with semidiameter of twice the Schwarzschild radius for the mass 
which has accumulated in the central black hole. 
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Fig. 12. Structure of the Paczynski-Wiita Model $10 in the equatorial plane near the end of the simulation. Panel a 
shows contours of the baryonic mass density (measured in gem -3 ) together with the flow field indicated by velocity 
vectors, panel b displays contours of the temperature (in MeV), panel c of the electron fraction Y e , and panel d of 
the entropy per nucleon. The length scale of the velocity vectors and the time elapsed since the beginning of the 
simulation of Model B64 of Ruffert & Janka (1998b) are given in the insert in the top right corner of the figures. 
The density contours are spaced logarithmically with intervals of 0.25 dex (bold: log(p [gem -3 ]) = 10.0, 11.0), the 
temperature contours are spaced linearly, starting with 1 MeV and 2 MeV and then continuing in steps of 2 MeV (bold: 
T = 10 MeV), the Y e contours are spaced linearly with intervals of 0.02 (bold: 0.02, 0.06, 0.10, 0.16), and the entropy 
levels were chosen to be 1, 2, 3, 4, 5, 6, 8, 10, 12, 14, 16, and 20 (bold: 6, 10, 20). Some of the contours are labeled with 
their respective values, and darker grey shading indicates higher values of a quantity. The circle around the center 
marks the inner "vacuum" boundary with semidiameter of twice the Schwarzschild radius for the mass which has 
accumulated in the central black hole. 
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Fig. 13. Contour plots of the mass fractions of a particles and heavy nuclei in the equatorial plane of Model B10 
(panels a and b, respectively) and of Model B10 (panels c and d, respectively) at the times given in the inserts in the 
top right corners of the figures. The contours are spaced logarithmically with steps of —1 (starting at —7; bold lines 
for levels —4 and —2) in case of the a mass fraction and with levels —8, —6, —4, —2, —1.5, —1, —0.8, —0.6, —0.4, —0.2 
and (bold lines for levels —6 and —1) in case of heavy nuclei. Typical mass numbers A of heavy nuclei are around 
80-120 with charge numbers Z w 25-50. Darker grey indicates higher values of the mass fractions. 



lows a lot of the gas which otherwise ends up in the torus 
because it has been spun off the surface of the rapidly 
rotating and oscillating merger remnant. The gas cloud 
that surrounds the compact central object has an average 
density above 10 12 gcm~ 3 initially, whereas at later times 
its density is below 10 12 gcm~ 3 because it was heated by 
shocks and viscous dissipation and becomes inflated due 
to the thermal gas pressure. If the simulation is started 



at a moment when this gas is in a post-merging phase of 
expansion, some of the matter (up to 2 x 10~ 2 Mq) can be 
ejected to leave the outer boundary of the largest grid and 
to become unbound (determined from the criterion that 
the total energy as the sum of kinetic, internal, and po- 
tential energies is positive). For the Newtonian Model B4 
which was carried on for a time long enough to follow this 
gas until it has reached the outer grid boundary, the value 
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Fig. 14. Azimuthally averaged radial structure of the accretion torus in the equatorial plane of Model B10 (left) and of 
Model B10 (right) at the end of the simulations. The curves show as functions of the equatorial radius d the logarithm 
of the density p(d) (measured in gem -3 ), the temperature T(d) (in MeV), the electron fraction Y e (d), the entropy per 
nucleon s(d), and the total cumulative mass in the torus, M(d) (in units of Mq). All curves are normalized to fit into 
the scale given on the ordinate and start at twice the Schwarzschild radius R s of the mass that has accumulated in 
the central black hole. Note that the gas density drops significantly inside d ss 3R S w 27 km. 



of the lost mass is in very good agreement with the mass 
l oss (0. 024 M ) found in Model B64 of Ruffert & Janka 
( 1998b ). This shows that the evolution of the gas swept 
out in large spiral arms becomes independent of the dy- 
namical evolution of the compact remnant at the center of 
the merger very early. In case of the Newtonian potential 
(left plot in Fig. ||), the thin solid lines labeled with M<j 
represent the mass on the grid with specific angular mo- 
mentum larger than the Kepler value at 3 Schwarzschild 



radii, 



3 



> 



3n 



v 



N 

Kepler 



{3R S )3R S = VGGM/c. For 



M the sum of the black hole and gas masses on the grid 
was used at all times. Since this overestimates the gravita- 
tional potential compared to the black hole in the numer- 
ical simulation, is too strict a limit on j. Therefore M$ 
is systematically somewhat smaller than M v but yields a 
reasonably good a priori estimate of the mass which can be 
found in the accretion torus in the quasi-stationary state. 

The maximum density on the grid as a function of 
time (Fig. [?]) shows a rapid drop from initially 2 to 
3 x 10 14 gcm~ 3 to about 30-100 times smaller values 
within only 1 ms, corresponding to the catastrophic dy- 
namical accretion after the start of the torus simulations. 
This steep decrease of the density is associated with the 
sudden replacement of the compact central object of the 
merger remnant by the vacuum sphere to model the black 
hole. The physical, more gradual generation of the event 
horizon will most probably soften this density drop some- 



what. The maximum density then levels off to a value 
which gradually decreases on the much longer evolution 
time associated with the subsequent accretion of torus ma- 
terial onto the black hole (Fig. |J). Since the mass remain- 
ing in the accretion torus is larger and its accretion time 
scale is longer in the Newtonian simulations, the value 
of the maximum density at the end of the simulations 
is higher there, and the temporal decrease of p max (t) is 
slower. The same tendencies can be seen for the maximum 
temperatures in the Newtonian and Paczyhski-Wiita mod- 
els (Fig. ||). The hottest spots in the accretion tori have 
temperatures around 10 MeV at the end of the simulations 
with a trend to somewhat larger values for the models 
where the black hole formation was assumed to occur late. 
In these cases the gas in the surroundings of the compact 
remnant of the merger has experienced additional heat- 
ing by the shocks and compression waves created by the 
violent oscillations of the central compact object. 

Since the temperatures are similar and the thermal 
energy of nondegenerate baryons contributes a dominant 
fraction of the internal energy of the gas, the specific in- 
ternal energies of the Newtonian and the Paczyhski-Wiita 
tori (plotted as functions of time in Fig. ||for Models B10 
and B10, respectively) are very similar, in both cases 
around 1.5 x 10 19 ergg _1 or about 15.5 MeV per nucleon. 
In contrast, the kinetic, internal, and potential energies of 
both models become very different towards the end of the 
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simulations (Fig. [T(]) which is explained by the large dif- 4. Properties of the accretion torus 

ference of the torus masses (compare M v in the two plots 

of Fig. pi). The structure of the accretion torus in the quasi 



stationary state is shown in Figs. 11 - 22[ In Figs. [IT] and 
|l2| contour plots of the density p, temperature T, electron 
fraction Y e , and entropy per nucleon s in the equatorial 
plane of the Newtonian Model BIO and of the Paczyhski- 
Wiita Model 2310, respectively, are given at the end of 
the simulations. Figure [l^ presents the corresponding in- 
formation for the mass fractions of a particles and heavy 
nuclei in both accretion tori. In Fig. [l4| the azimuthally av- 
eraged radial structure of the tori in the equatorial plane 
(quantities p, T, Y ei and s) is displayed, together with the 
cumulative torus mass M{d) as function of the equatorial 
radius d. Contour plots for p, T, Y e , and s in the x-z- and 
y-z-planes perpendicular to the equatorial plane are shown 
in Figs. p"5| and |l6| . Integral mass distributions perpendic- 
ular to the equatorial plane for the initial model of our 
simulations as well as for the final states of Models BIO 
and 2310 are presented in Fig. |l7|. Finally, Figs. 18-22] give 
information about the azimuthal and radial velocities and 
the specific angular momentum of the gas in the equatorial 
plane of the initial model and of the evolved tori. 

At the end of the simulations the tori have become 
nearly axially symmetric with only minor deviations (Figs. 
|TT| and |l2]). Two hot spots can be seen close to the in- 
ner grid boundary at 2R S which coincide with density 
maxima. They continue to carry the memory of the two 
very prominent spiral-like arms which are formed dur- 
ing the neutron star merging, grow right afterwards, and 
are wound up into the toroid around the black hole dur- 
ing the subsequent evolution. The torus of Model 2310 in 
Fig. |l2| is smaller — the p — 10 10 g cm~ 3 contour is at 70- 
80km — than the torus of Model BIO in Fig. [fl] where 
the p = 10 10 g cm -3 contour extends out to 120-130 km. 
There are two reasons for that. On the one hand, the gas 
mass which remains around the black hole is smaller in 
the former model (see Table [l]), on the other hand the 
gravitational potential is stronger in the Paczyhski-Wiita 
case. This is the reason why despite of a difference of a 
factor 7 in the torus mass, both models have average den- 
sities which differ only by a factor 3 (about 3 x 10 11 g cm~ 3 
for Model B10 compared to approximately 10 11 gcm~ 3 in 
case of Model 2310, see Fig. [l4]). The entropies in both 
models are nearly the same (between 5 and 10 fee per nu- 
cleon) and the entropy profiles are very similar in the re- 
gion where most of the mass is sitting (Fig. fu]). The New- 
tonian torus is somewhat hotter, its average temperature 
is around 6MeV compared to 4MeV for Model 2310. It is 
also more neutron rich with a mean value of Y e w 0.05- 
0.15. Model 2310 has typical Y e values between 0.1 and 
0.17 because of its lower density. This has the consequence 
that the neutrino opacities are smaller and thus the in- 
crease of the proton abundance by the emission of elec- 
tron antineutrinos (and the cooling by the loss of neutri- 
nos and antineutrinos of all flavors) proceeds faster. Since 
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the hot, expanded neutron star matter radiates predom- 
inantly electron antineutrinos (see Sect. ^J), the initially 
very neutron rich state evolves towards higher electron 
and proton number fractions. Due to the high tempera- 
tures, electrons are not degenerate in the tori and therefore 
positrons are abundant, in considerable numbers. Typical 
electron degeneracy parameters ry e = /i e /T (/i e is the elec- 
tron chemical potential) are around 2, some regions have 
values of about 4, and in large regions one finds rj e < 2. 
For the conditions in the tori, in particular due to the 
rather high entropies, the nucleons are mostly unbound 
and a particles and heavy nuclei are present only in small 
numbers. The maximum mass fractions X a of a particles 
are around a few per cent, and heavy nuclei appear in 
significant abundances only where the temperature drops 
below IMeV (Fig. [§. 

The perpendicular cuts (Figs. |l5| and |l6|) confirm the 
nearly axially symmetric structure of the tori in Mod- 
els BIO and BIO at the end of the simulations. Again, 
some differences between the x-z- and y-z-cuts reflect the 
last remainders of the spiral arms which have been in- 
flated and dissolved into the tori. While the temperature 
and density contours show a rather regular shape, primar- 
ily determined by the balance of pressure gradients and 
gravitational and centrifugal forces, the electron fraction 
and entropy are more irregular and patchy because both 
quantities carry information about the whole preceding 
evolution, in particular about the integral effects of neu- 
trino emission and non-adiabatic hydrodynamic processes. 
Near the poles of the black hole and along the system axis, 
the density has decreased to values below 5 x 10 8 gcm~ 3 . 
This is only one order of magnitude above the lower den- 
sity limit which is set to 5 x 10 7 g cm -3 in the surround- 
ings of the torus for numerical reasons, but nevertheless it 
is more than 3-4 orders of magnitude below the average 
densities inside the tori. In this sense we see the formation 
of an "evacuated" , cylindrical funnel along the rotational 
axis of the black hole-torus system. Material which was 
swept into the polar regions during and immediately af- 
ter the merging of the neutron stars falls into the newly 
formed black hole very quickly within a free-fall time scale 
because it is not supported by centrifugal forces. A com- 
parison of the initial condition for the torus simulations 
of Models BIO and BIO (corresponding to the si tuation 
at t « 10 ms in Model B64 of Ruffert & Janka |l998b| ) 
with the quasi-stationary states about 5 ms later reveals 
this rapid cleaning of the axial funnel (see Fig. [Tt]). In the 
upper panels of the three figures in Fig. [l?] the contours 
contain all points (d, z) where the cumulative gas mass 
M gas(d,20 given by 



M'{d,z) 



d+Ad 



2tt 



AUd / ckpp(r) 



is constant (the labels at the contours represent logarith- 
mic values of the mass measured in Mq). The integration 
is done on a cylindrical grid with coordinates (d, z, tp). The 



integral of Eq. (|3|) therefore sums up all the mass within a 
hollow cylinder of thickness Ad = 1 km extending from z 
to infinity (in praxi: the upper grid boundary). The con- 
tours in Fig. [H] are mirrored along the system axis at 
d = 0. Note that only the mass on one side of the equato- 
rial plane is added up. The lower panels in Fig. ^ show 
the contours that correspond to constant values M gas (d, z) 
according to the integral 



M gas (d,z) 



2tt 



(4) 



lo Jo 

which sums up the gas mass in a cylinder with radius d 
that extends from z to infinity around the system axis. The 
panels in Fig. ^ give detailed information about the mass 
distribution in the surroundings of the accretion torus. Ini- 
tially, about 10~ 3 Mq of gas were distributed above the 
compact massive object at the center of the merger (se e 
second panel from the top in the left column of Fig. |l7| ). 
However, with increasing vertical distance z from the 
equatorial plane the mass drops extremely rapidly already 
in the post-merging configuration. In the final states, the 
total mass inside a cylinder with radius d ps 15 km is only 
a few 10 -4 Mq, most of this gas is very close to the equa- 
torial plane. Around the rotation axis the quasi-stationary 
states of Models BIO and £>10 look very similar. The larger 
torus mass of the Newtonian computation, however, leads 
to a higher mass concentration near the equatorial plane 
at radii d ^ 20 km. Differences are therefore visible in this 
region for the contours corresponding to cumulative gas 
masses above 10 -5 Mq. We note here that the simulations 
described in this paper do not include the effects of neu- 
trino energy deposition in the vicinity of the black hole. 
For this reason and because we do not allow the gas den- 
sities to drop below 5 x 10 7 g cm -3 on the grid, our models 
most likely overestimate the mass of the gas surrounding 
the torus. 

Information about the motion of the gas in the initial 
and final models of the torus simulations is provided by 
Figs. [l£ - 22 . In Fig. |l8| the azimuthal velocities v v (d) are 
plotted for all grid zones in the equatorial plane versus the 
distance d from the grid center (dots) immediately before 
the black hole is assumed to form. The spread of the points 
at a given radius d reflects the deviations from rotational 
symmetry. For an axially symmetric configuration all dots 
at a specific radial distance would cluster on top of each 
other. Also, the motion of the gas is clearly sub-Keplerian 
which indicates the importance of pressure support in the 
object. This can be seen from a comparison with the bold 
solid line which represents the local Newtonian Kepler ve- 
locity, 



(3) ^epler(rf) 



GM (d) 



(5) 



given as a function of the equatorial radius d, with M(d) 
being the mass enclosed by the sphere of radius d. The sit- 
uation is different at the end of the computed torus evolu- 
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tion about 5 ms after the assumed formation of the black 
hole (Fig. [II] for the Newtonian simulation and Fig. ^ for 
the Paczyhski-Wiita case) . The spread of the dots has de- 
creased, indicating that the torus is much more isotropic 
in ip than the post-merging configuration. For distances 
d ^ 30 km the Newtonian torus has azimuthal velocities 
larger than the local Keplerian value. These allow the gas 
between 2 R s and 3 R s to remain on orbits around the 
black hole despite of a positive density and pressure gra- 
dient in this region (see Fig. |l4|). Also for d > 30 km the 
pressure support is important. This is suggested by the 
significant drop of the azimuthal velocities below the Ke- 
plerian values in Fig. |l9|. In the torus of Model BIO the 
orbital velocities are much closer to the Kepler velocity 
for the Paczyhski-Wiita potential. The value of 



of electron neutrinos (f e ), electron antineutrinos {y e ) and 



v 



PW 
Kepler 



' GM(d)d 
(d-R s ) 2 



(6) 



is twice as large at d = 2R S than its Newtonian counter- 
part, ^K^ ler (2i? s ) = 2wK eplor (2i? s ) = c, and therefore the 
orbital velocities within d rs 50 km are significantly larger 
in Model BIO than in Model BIO. 

The specific angular momentum j(d) = v ip {d)d of the 
matter in the equatorial planes of both models has typi- 
cal values of (3-6) x 10 16 cm 2 s _1 (Fig. |l|). In the New- 
tonian torus the lines for j(d) and the Keplerian value 
JKe P ier( d ) = w Ke P ior ( d )d intersect at d » 30 km rs 3ii s 
which is roughly the position of the density and pres- 
sure maximum (Fig. |l4|). Inside this radius one has j{d) > 
ixepierC^)- ^ n contrast, for the Paczyhski-Wiita potential 
the curves of j(d) and j™ p x eT {d) = v™ plcl (d)d touch at 
d ~ 3i? s where iKepler(^) uas a minimum corresponding 
to the last stable circular orbit. 



Figure 22 shows the radial velocities for all grid zones 



in the equatorial plane. The Newtonian Model BIO has 
achieved a quasi-stationary state with very small radial 
velocities v r at the end of the simulation. The average 
inflow velocity between 25 km and about 130 km is (v r ) rs 
— 10 8 cms -1 . At smaller distances from the black hole the 
gas is rapidly falling in, at distances beyond 130 km the 
dilute outer parts of the disk are slowly expanding due 
to the outward transport of angular momentum in the 
torus. In contrast, the Paczyhski-Wiita Model £510 is still 
evolving and has not developed stationary conditions. It 
expands for 40 km ^ d Si 80 km with radial velocities up 
to 5% of the speed of light whereas the gas interior to d rs 
3i? s rs 27 km is collapsing very rapidly into the black hole, 
and also the dilute gas exterior to 80 km moves inward 
with large velocities. 

5. Neutrino emission 

The neutrino luminosities as functions of time for the New- 
tonian and Paczyhski-Wiita tori are displayed in Figs. ^3] 
and E4L respectively, and the corresponding mean energies 



heavy-lepton neutrinos (y x 



v T , v T ) are given in 



Figs. |25| and |26j. In order to show the sensitivity of the 
results to the moment the black hole is assumed to form, 
data from two simulations which were started at different 
times are plotted in each of the figures: Models B4 and 
B10 for the Newtonian simulations, and Models B2 and 
£>10 for the Paczyhski-Wiita potential. 

Despite of its much smaller mass (by a factor 7.5) 
and lower density and temperature (see Fig. |l4|), the 
Paczyhski-Wiita torus radiates neutrinos with similar lu- 
minosities and mean energies as the Newtonian model. 
Mainly v e and v e are emitted. At the end of the computed 
evolution (t rs 15ms), the Newtonian model (Figs. 23 
p5| ) has a total luminosity L v = 12 x 10 52 ergs 
contributions of L„ e = 3.5 x 10 52 ergs _1 from v P , L 



and 
with 



6.5 x 10 52 ergs _1 from 9 e , and L Vx = 0.4 x 10 52 ergs _1 
from v x individually. The mean energies of the emit- 
ted neutrinos are (e^) = 9MeV, (ep e ) = 13MeV and 
(e^) = 21MeV. For the Paczyhski-Wiita torus, the total 
luminosity reaches about 50% of the value obtained in the 
Newtonian simulation, L v — 6.7 x 10 52 ergs _1 , and the 
luminosities of v e and v e reach 60-70% of the correspond- 
ing Newtonian values. The average energies of v e and v e 
are slightly higher (around 10-12 MeV and 15-17 MeV, 
respectively). The differences between the simulations are 
more pronounced in case of v x . At t rs 15 ms the v x lu- 
minosity is approximately one order of magnitude smaller 
for the Paczyhski-Wiita model, and the emitted v x are less 
energetic with a mean energy of only 16-19 MeV instead 
of 21 MeV. 

These differences result from the fact that the 
Paczyhski-Wiita torus is essentially transparent for muon 
and tau neutrinos whereas a well defined average muon 
and tau neutrinosphere (to be more precise: a toriodal 
neutrinosurface) exists in the Newtonian model. This is 
clearly visible from the dotted lines in panels c of Figs. [27] 
and |28|, which represent the intersections of the equato- 
rial plane at z — with the neutrinosurface. The latter 
is defined as the two-dimensional hypersurface where the 
optical depth of the torus perpendicular to the equatorial 
plane is unity, i.e., where 



z , Vi {x,y) = / dz' K Vi (x,y, z') = 1 , 



(7) 



with K Ui (x,y, z) being the total opacity (defined as the 
inverse of the mean free path) for the energy transport 
of neutrino v± at a point (x, y, z). In contrast, the very 
similar properties of the v e and v e emission can be under- 
stood from similar thermodynamical conditions (density, 
temperature, entropy, see Fig. [HJ) at the corresponding 
neutrinosurfaccs. The latter have an outer radius of about 
100 km in case of the Newtonian torus (panels a and b in 
Fig. |27]) and of 70-80 km in the less massive Paczyhski- 
Wiita model (panels a and b in Fig. ^8|). The slightly dif- 
ferent sizes and thus different areas of the toroidal neutri- 
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nosurfaces in panels a and b of Figs. ^ and ^ account 
for the moderate differences of the v e and v e luminosities 
in both models. 

The neutrino emission is primarily determined by the 
size of the neutrinosurface and the thermodynamical con- 
ditions in the layer where the average neutrino optical 
depth is around unity. For this reason the torus mass influ- 
ences the luminosities indirectly through the radius of the 
neutrinosurface, until the torus mass and its density and 
temperature become so low that neutrino transparency is 
reached. The smaller torus mass and larger relative and 
absolute accretion rate in the Paczyhski-Wiita simulation 
(see Figs. |, | and |) lead to a different time evolution 
of the neutrino luminosities (and mean neutrino energies) 
in both models. Extrapolation of the luminosity decrease 
towards the end of the simulations in Figs. ^ and sug- 
gests a longer decay time scale for the Newtonian model 
in agreement with the longer accretion time scale given in 
Fig. |. 

Comparing the models with early formation of the 
black hole, Model B4 and Model B2, with those where 
the black hole collapse is assumed to happen later, Mod- 
els BIO and BIO, shows that the latter have somewhat 
higher neutrino luminosities. This is explained by the 
higher temperatures of the tori of Models BIO and £510, 
see Fig. ||, and the correspondingly lower densities in 
the thermally inflated later states (Fig. (7J). The effect is 
particularly strong for v x which are mainly produced by 
the annihilation of electrons and positrons into neutrino- 
antineutrino pairs at the conditions present in the tori, 
because the energy emission rate for this process is ex- 
tremely temperature sensitive and increases proportional 
to T 9 . Since the accretion rates are time dependent and 
the temperatures in the tori show fluctuations, the neu- 
trino luminosities and mean energies are variable on a time 
scale of 1-2 ms. 

The regions with the strongest neutrino emission are 
visible as dark grey shaded areas in 



and 



In 



the former two figures the contours correspond to levels 
of constant energy emission rate per unit area as obtained 
by integration of the energy loss rate per volume from the 
equatorial plane at z = to infinity. In Fig. |2^ vertical 
cuts through the Newtonian and Paczyhski-Wiita accre- 
tion tori are displayed which show the total neutrino en- 
ergy loss rates per unit volume in the x-z and y-z planes. 
The peak values of the emission rates are similar in both 
types of models but the main neutrino emitting region 
is less extended in case of the Paczyhski-Wiita potential. 
The dotted lines in Figs. and mark the intersection 
of the equatorial plane with the two-dimensional hypersur- 
face where the transport optical depth for the energy flux 
is unity (Eq. ^), i.e., where the spectrally averaged mean 
free path of the neutrinos or antineutrinos of a certain fla- 
vor is of the same order as the size of the emitting volume. 
Outside this neutrinosurface the neutrinos stream off es- 
sentially freely and interact by scattering or absorption 



with the gas particles on average only one more time. The 
muon and tau neutrino opacity is dominated by neutral- 
current scatterings off neutrons and protons. Only v e and 
v e are also absorbed on neutrons and protons, respectively, 
via charged-current inverse beta processes. The neutri- 
nosurface of the heavey-lepton neutrinos has a toroidal 
shape only in case of the Newtonian calculation (Fig. 
panel c) but splits up into several distinct islands for the 
less massive and less dense Paczyhski-Wiita torus which is 
near to neutrino-transparent conditions. The highest neu- 
trino energy loss rates are found in a region between the 
inner grid boundary at 2R S w 18 km and an equatorial 
radius of approximately 70 km. In the Newtonian model 
significant contributions to the neutrino luminosity come 
even from larger distances out to about 100 km where the 
larger volume compensates for the smaller emission rates. 
The vertical cuts of Fig. |2^ confirm an effect which was 
already visible in the density plots of Figs. [lB] and |l6|: In 
the Paczyhski-Wiita model, in contrast to the Newtonian 
simulation, gas flows towards the black hole even from 
the poles where it shows up by its neutrino emission in 
the right plot of Fig. |2^. 

The maximum energy loss rates of neutrinos and an- 
tineutrinos of all flavors are typically of the order 10 32 erg 



,-3 -i 



,-3 -i 



In single peaks values of even 10 erg cm s 
can be reached (Fig. |2{J). In gas with density p 



10 11 gcm 3 (Figs. |15| and |16|) these rates correspond to 



a specific energy loss of 1000 MeVs -1 per nucleon up to 
even 10 4 MeV s _1 per nucleon just before the gas reaches 
the inner grid boundary and disappears in the black hole. 
For a total neutrino luminosity of 10 53 erg s _1 from a torus 
with mass of approximately 0.25 Mq in the Newtonian 
simulation, one calculates an average neutrino energy loss 
rate of 200 MeVs" 1 per nucleon. This means that the 
binding energy of a nucleon in the gravitational potential 
of the 3 Mq black hole at the position of the last stable 
circular orbit at 3i? s ~ 27 km is radiated away in less than 
a second, or an energy equivalent of more than 2% of the 
nucleon's rest mass escapes in neutrinos within only 0.1 
seconds. This estimate is in agreement with the instan- 
taneous efficiency for the conversion of rest-mass energy 
into neutrino energy, q v = L v / '(Mc 2 ), which we calcu- 
late at the end of the simulation to be « 1.3% for 
the Newtonian torus (Table |J) . In case of the Paczyhski- 
Wiita model the corresponding value is about a factor of 
2.5 lower, g^ w w 0.5%. These numbers are significantly 
smaller than the maximum efficiency of 5.7% for relativis- 
tic disk accretion onto a Schwarzschild black hole (8.3% 
for a thin Newtonian accretion disk) because the tori are 
not transparent for neutrinos. Due to the high densities 
and temperatures, the diffusion time scale for neutrinos 
becomes longer than the accretion time scale of the gas 
into the black hole. Therefore the tori are advection dom- 
inated and cooling does not reach its maximum possible 
efficiency. 
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Table 2. Torus mass M v , black hole mass accretion rate M v , typical accretion time scale of the torus, £ acc = M v /M v , 
total neutrino luminosity L„, efficiency q v = L v / '(Mc 2 ), integral rate of energy deposition by neutrino-antineutrino 
annihilation around the accretion torus, E vi/1 and efficiency = E v a / L v for the Newtonian Model BIO and Paczyhski- 
Wiita Model 2510. All quantities are given at the time ifi n when the simulation was stopped. E vD gives the estimated 
total energy deposition by vv annihilation in the time interval tg n + i a cc which is roughly equal to the duration of 
neutrino emission. 
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6. Neutrino-antineutrino annihilation 

We evaluate our hydrodynamical models in a post-process- 
ing step for the annihilation of neutrinos and antineutrinos 
into electron-positron pairs. The correspon ding m ethods 
were explaine d in d etail in Ruffert & Janka ( 1998a ) and in 
Ruffert et al. ( 1997 ). Neutrinos and antineutrinos emitted 
from the hot accretion torus interact with each other in 
the surroundings with a finite probability which depends 
on the number densities and energies of these neutrinos 
and on the angle between the directions of neutrino and 
antineutrino propagation (see Goodman et al. 1987, Coop- 
erstein et al. |1987| , Mayle |199C| ; also Ruffert et al. |1997| ). 
Therefore the local energy deposition rate by vv annihila- 
tion increases proportional to the product of neutrino and 
antineutrino luminosities and the spectrally averaged neu- 
trino energy, times a factor that accounts for the depen- 
dence on the angular distribution of the neutrinos. The 
annihilation rate drops rapidly with increasing distance 
from the neutrino source because both the neutrino num- 
ber densities and the mean angle of neutrino-antineutrino 
collisions decrease. 

The computational procedure to obtain the energy de- 
position outside the accretion torus can be briefly sum- 
marized as follows. In a first step the hydrodynami- 
cal model at a chosen time is mapped from the nested 
grids of the simulation onto an equidistant cartesian grid. 
Next, for each type of neutrino or antineutrino the two- 
dimensional surface is determined where the optical depth 
in z-dircction is unity (see Eq. |^). The neutrino energy 
loss rates from the torus volume are projected onto these 
surfaces and treated as surface emissivities. In order to 
compute the energy deposition rate per unit volume by 
vv annihilation at a point r, one has to sum up the con- 
tributions by the neutrino and antineutrino emission from 
all parts of the neutrinosurfaces which radiate into the 
direction of point r. The total energy deposition rate in- 
cludes contributions from all three flavors of neutrinos and 
antineutrinos. Finally, volume integrals of the energy de- 
position rate can be obtained by summation over specified 
regions of the equidistant cartesian grid. The results of this 
evaluation are plotted in Figs. |3(] and [H] for the final states 



of our torus simulations. The figures show azimuthally av- 
eraged (around the z-axis) quantities in a plane perpen- 
dicular to the equatorial plane. 

Figure ^ displays contours of constant energy depo- 
sition rate per unit volume in those regions around the 
Newtonian (left) and Paczyhski-Wiita (right) tori where 
the baryonic mass density is less than 10 11 gcm~ 3 . The 
dotted lines represent density contours, the dashed lines 
mark the positions of the average neutrinosurfaces of v e , 
v e and v x (from outside outward) . The v e neutrinosurface 
is very close to the toroidal surface which corresponds to a 
density of p = 10 11 g cm~ 3 . The neutrinosurfaces of v e and 
v x are deeper inside, because v e are absorbed onto protons 
which are less abundant in the torus than neutrons, and 
the streaming of the v x is inhibited only by neutral-current 
scatterings off nucleons. Both the dashed and dotted sets 
of contours demonstrate that the Paczyhski-Wiita torus 
is significantly smaller and nearly transparent for electron 
antineutrinos and heavy-lepton neutrinos. 

The maximum energy deposition rates by vv annihila- 
tion (solid contours in Figure^) exceed 10 30 ergcm^ 3 s~ 1 
in the polar regions above and below the equatorial plane 
at heights \z\ between 10km and 30km in the Newtonian 
model. Such high values occur also in the Paczyhski-Wiita 
simulation, but only very close to the surface of the torus 
and to the equatorial plane, between the black hole and 
the last stable circular orbit at 3R S ~ 27 km. the values 
at larger distances \z\ are typically one order of magni- 
tude lower than in the Newtonian model at the same \z\. 
Within the displayed region, the energy deposition rate 
along the polar axis decreases roughly proportional to z~ 6 
for 1 2: | > 30 km. 

The integral rate of heating by vv annihilation in the 
computed volume is 4.9 x 10 50 ergs _1 for the Newtonian 
model. About 40% of this energy (1.9 x 10 50 ergs _1 ) are 
deposited in a cylinder with radius 3i? s = 27 km around 
the z-axis, 90% (4.4 x 10 50 ergs _1 ) in a cone with open- 
ing half-angle of 30-40 degrees which is approximately 
bounded by the isodensity contour for p = 10 10 gcm -3 
in Fig. |3^. The corresponding numbers for the Paczyhski- 
Wiita torus are 3.1 x 10 50 ergs^ 1 for the total energy de- 
position, of which 1.3 x 10 50 ergs _1 end up in the cylin- 
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Fig. 31. Cumulative energy deposition rates by vv anni- 
hilation for the Newtonian Model B10 at time t = 14.9 ms 
(upper panel) and for the Paczyhski-Wiita Model $10 at 
time t = 15.2 ms (lower panel). The azimuthally averaged 
annihilation rate per unit volume was integrated along the 
z-direction from a given value of z out to the grid bound- 
ary (practically infinity because of the rapid decrease of 
the annihilation rate with distance) and in addition from 
d = to the value of d given on the abscissa. The contours 
are spaced logarithmically in steps of 0.5 dex and repre- 
sent values only for the hemisphere above the equatorial 
plane (measured in ergs^ 1 ). They are mirrored along the 
z-axis. The two panels should be compared with the lower 
panels in Fig. [l^ and were constructed in an analogous 
way. 

der along the axis and 1.9 x 10 50 ergs _1 in the cone. 
From these values one computes the total conversion ef- 
ficiencies q V i, = EyjjjLy of neutrino energy into energy 

of the pair-plasma fireball to be q^ fJ = 4.1 x 10~ 3 and 
g PW = 4 6 x 1(r 3 ( Tab i e |). 

Figure |3l] provides information about the spatial dis- 
tribution of the cumulative energy deposition rates in both 
models. The plots show z-d-integrals of the rates per unit 
volume, <5i/p(r), which are computed as 



E vl/ (d, z) 



d X / d££ / d^Q„ p (r) . 
Jo Jo 



(8) 



These plots are thus constructed in analogy to the panels 
of the z-d integrals for the mass in Fig. [l?]. Note that the 
values correspond to the energy that is deposited only on 
one side of the equatorial plane. A large part of the en- 
ergy, nearly 10 50 ergs _1 , is dumped within 30 km around 
the polar axis and at heights larger than about 30 km 
above (and below) the equatorial plane in the Newto- 



nian model and at heights larger than about 10 km in the 
Paczyhski-Wiita case. According to Fig. [l7] there are less 
than 10" 4 M Q of gas in these regions in both models. A 
comparison of the two panels in Fig. ^ reveals that in 
the Paczyhski-Wiita simulation the energy is more con- 
centrated towards the equatorial plane so that at heights 
\z\ 30 km the numbers are about one order of magnitude 
lower than in the Newtonian model. This difference is ex- 
plained mainly by the smaller size of the neutrinospheres 
in the less massive Paczyhski-Wiita torus. A minor part 
of the effect is also due to the lower neutrino luminosities 
which reach only little more than half the values of the 
Newtonian model. 



7. Analytical estimates 

7.1. Viscosity and evolution 

The evolution of the Newtonian torus on longer time scales 
will be governed by the outward-directed viscous transport 
of angular momentum. Since the Euler equations which 
are solved numerically with the PPM method do not con- 
tain viscosity terms, and since this scheme does not require 
any artificial viscosity to treat shock waves, the most im- 
portant viscous dissipation in the absence of shocks comes 
from the numerical viscosity of the code. The latter is as- 
sociated with the discretization of the equations and thus 
depends on the chosen grid resolution, but is also deter- 
mined by the input physics implemented in the hydrody- 
namics code. In case of the quasi-stationary Newtonian 
Model B10 we shall estimate the size of this numerical 
viscosity from the torus properties. 

With the average value of the radial inflow velocity, 
(v r ) w — 10 8 cms^ 1 (see Fig. |2| and Sect. ||), and a mean 
radius (d) of about 50 km we estimate an accretion time 
scale of t acc = (d)/(v r ) rj 50 ms in very good agreement 
with the independent determination via the mass accre- 
tion rate of the black hole, t acc = M v /M v w 53 ms (see 
Fig. |5| and Table ||). The viscous time scale is given by 
t vis w M v /(6TrnR s ) — see Eq. (13) in Ruffert et al. fll997| ) 
- with the dynamic viscosity n ~ a(pH URepier) where 
H ~ R s is the height of the torus, p the average den- 
sity in the torus, WKcpior = ^Kc P \cM) d = y/GM^a/d the 
Keplerian velocity at a representative radius d ~ 4i? s , 
and a the dimensionless a-viscosity parameter. Using 
p = M v /Vt with the torus volume Vt ~ 2ir 2 R^d, one finds 
a ~ [SlKc P icr(d)£vi s r ■ Setting t vis = i acc = (d)/(v r ) we 
get a value for the a-parameter associated with the nu- 
merical viscosity of a ~ (v r ) /fKepier(rf) ~ 0.01. It is inter- 
esting to note that the torus shapes computed by Popham 
& Gammie ( 1998 ) for such values of the a- viscosity look 
very similar to the cross sections through the tori of our 
simulations plotted in Figs. [l5| [l6] and |3^. 

The numerical viscosity in our simulations therefore 
corresponds to a value a ~ 0.01, which is a bit higher than 
the "optimum" value a* — rf / (p-fffKepier) where the vis- 
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cous energy dissipation and the energy emission by neutri- 
nos are balanced, a requirement which ensures maximum 
efficiency for the conversion of rest-mass energy into neu- 
trinos. Ruffert et al. (|1997|) estimated a* ~ (1.7... 6.6) x 
10" 3 (i? s /9km) 7 / 2 (M v /0.1 M )" 3 / 2 on grounds of a very 
simple one-zone model of the torus. Consistent with the 
larger value of the viscosity, the mass accretion rate of 
the black hole was found to be M v « 5M Q s _1 in the 
numerical simulation (Fig. || and Table ^), which is some- 
what higher — and the corresponding lifetime < acc of the 
torus somewhat shorter — than estimated by Ruffert et 
al. (1997) for the optimum value a* . This implies that the 
torus does not lose energy in neutrinos at the maximum 
theoretical efficiency of about 8.3% for disk accretion on 
a nonrotating black hole in Newtonian gravity. I.e., the 
torus in the numerical simulation is advection dominated, 
in agreement with the findings in Sect. (see also Ta- 
ble g). Because of the high densities and temperatures in 
the torus, the matter is not transparent to neutrinos and 
the neutrino diffusion time is longer than the accretion 
time of the gas. Therefore a sizable fraction of the gravita- 
tional binding energy that is dissipated into heat is carried 
into the black hole before neutrinos are able to transport 
it away. Although the neutrino luminosities obtained in 
the numerical simulations are close to those estimated by 
Ruffert et al. (1997), the shorter lifetime of the torus leads 
to a smaller time integral of the radiated energy. 

Popham et al. ( 1998 ) have calculated self-consistent 
stationary models for neutrino emitting accretion disks of 
black holes for a wide range of parameters. With the large 
accretion rates and high densities, the tori obtained in 
our simulations are in an extreme corner of their param- 
eter space where the idealized treatment of the neutrino 
cooling by Popham et al. ( 1998[ ), who essentially assumed 
neutrino-transparent conditions, reaches its limits of ap- 
plicability. 

Finally, we point out that neutrino emission is rather 
inefficient in transporting away angular momentum and 
therefore this is not the driving force of the torus evolution 
in our investigated cases. Setting the mass accretion rate 
by the black hole in relation to the angular momentum loss 
rate, J ~ M v fl^ cp i OY d 2 , and requiring J to be equal to the 
angular momentum carried off by neutrinos, J = J v ~ 
[L v I c 2 )£lKepierd 2 , one finds M v ~ L v jc 2 . For a neutrino 
luminosity of L v — 10 53 
rate of the order of M v ~ ^qMq s 
of magnitude below our numerically determined numbers. 



erg s 1 this gives a mass accretion 
1 which is two orders 



7.2. General relativistic effects 

The torus simulations described in this paper were per- 
formed with basically Newtonian physics except that the 
gravitational potential was replaced by a Paczyhski-Wiita 
potential in some cases. Also the emission of gravitational 
waves and the corresponding back-reaction on the hydro- 
dynamic flow were included, which, however, is not im- 



Table 3. In the second and third column ratios of gen- 
eral relativistic ("GR") to Newtonian ("N") results are 
given for the integral energy deposited by vv annihilation 
for the different values of a = Rt/R s listed in the first 
column. The other columns contain values of the factor 
U/(af3V) in Eq. (pl|). This factor decides about the im- 
portance of neutrino-electron and neutrino-positron scat- 
tering relative to vv annihilation. It is assumed that the 
neutrinos are emitted from a toroidal neutrinosurface with 
radius Rt and center at 3R S + R t . The energy deposition 
rates by annihilation or scattering are integrated either 
over the volume of a cylinder with radius 3i? s around the 
system axis, or over the volume of a cone enclosed by the 
rotated boundary curve of Eq. (flO|), respectively. 
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portant for the nearly axially symmetric black hole-torus 
configurations. 

The influence of general relativistic effects on the neu- 
trino emission and on the neutrino-antineutrino annihi- 
lation are not obvious because of several competing ef- 
fects which partly cancel each other (see, e.g., Shapiro & 
Teukolsky 



1983, Luminet 197S, Marck 1996, Cardall & 



199S). The effects considered here 



Fuller |1997| , Usui et al. 
are : 

(1) General relativistic blueshifting of those neutrinos 
which are emitted from the accretion torus in the direction 
of the black hole before they annihilate, and redshifting 
of those neutrinos which annihilate at distances from the 
black hole larger than the radius of their production. 

(2) Gravitational redshifting of the energy as the electron- 
positron pair plasma expands away from the site of its 
creation by vv annihilation. 

(3) General relativistic ray bending of the neutrino trajec- 
tories which changes the average angles at which neutrinos 
and antineutrinos collide and thus also changes the rate 
of vv annihilation. 

(4) Gravitational trapping of the energy which is deposited 
in the close vicinity of the event horizon by vv annihila- 
tion and therefore is not able to escape from the strong 
gravitational attraction of the black hole but is ultimately 
sucked in by the hole. 

It is not easy to implement these effects consistently 
into the hydrodynamic simulations and into the post-proc- 
essing procedure which we used to evaluate our models 
for neutrino-antineutrino annihilation. In order to esti- 
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mate the relative changes due to general relativity ap- 
proximately, we therefore consider a simplified treatment 
for an idealized model geometry which mimics sufficiently 
closely the situation described by the hydrodynamic re- 
sults. For this purpose, we assumed that the neutrinos are 
emitted from a ncutrinosurface that has the shape of a 
torus centered on the equatorial plane. At least in case of 
the Newtonian accretion disk, which is optically thick for 
neutrinos, Fig. |3^ shows that this picture is indeed appro- 
priate. In the following analysis we shall explicitly presume 
that the black hole-torus configuration is axially symmet- 
ric. We neglect effects due to the relativistic Doppler shift 
of the neutrinos emitted from the torus although the gas 
orbits around the black hole at a fair fraction of the speed 
of light. 

Integrating the energy deposition rate by vv annihila- 
tion over a volume V outside of the torus, one can write 



d 3 rQ^(r) = R- 1 -KCV, 



(9) 



where TZ = 9RiR^ 2 (3R s + R t )~ 2 is a dimensionless fac- 
tor that depends on the torus radius R t , the horizon ra- 
dius of the black hole, R s , and the radius of the torus 
center, 3i? s + Rt- The factor C contains the weak interac- 
tion coefficients c„ iPi (i = e, /i, r) and is proportional to 
the product of the neutrino and antineutrino luminosities 
multiplied by the sum of the average neutrino and antineu- 
trino energies which determine the vv annihilation rate, 
i.e., C oc J2i c v i v i L Vi L i j i ((e Vi ) + (ep 4 ))- The term V is also a 
dimensionless factor that results from the volume integral 
and carries information about the angular distributions of 
the neutrinos and antineutrinos at the positions r. 

Two different axially symmetric volumes are consid- 
ered. On the one hand, the integration is performed over 
a cylinder with radius 3R S around the symmetry axis of 
the torus-black hole system. On the other hand, the vol- 
ume of a cone is considered whose surface is described by 
a rotation of the curve 



> i (10) 




for > 1 

3R. 



around the symmetry axis. In all points on this surface the 
component of the gravitational force parallel to the equa- 
torial plane is balanced by the centrifugal force of mat- 
ter which has angular momentum equal to the Keplerian 
value at radius 3R S , j = jKepier (3Rs) = VQGMbh/c. This 
surface may approximate the boundary of a baryon-free 
funnel along the axis. 

If one assumes that the value of the integrand on the 
axis at vertical height z is representative of all values in 
the plane defined by this z, one can make use of the axial 
symmetry which allows one to calculate the vv annihila- 
tion rate along the axis particularly easily. The integra- 
tions parallel and vertical to the system axis can then be 



separated, and only a one-dimensional integral for V is left 
for numerical evaluation : 



(11) 



First we consider the cases without general relativis- 
tic effects. For the cylindrical volume one sets £ = z/R s 
and starts the integration at z = R s , i.e. £ m j n = 1. The 
integrand is found to be 
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(13) 



where a = Rt/R s and f3 = 3 + a. The minus sign in 
Eq. ( p"3| ) corresponds to fiQ — cos 9q , the plus sign to fii = 
cos 6*i, if 6q > 6i and 9 and 9\ are the angles between 
the symmetry axis and the direction of rays that leave the 
torus surface tangentially and go to an observer on the 
axis. Taking the variable of integration to be £ = xj (3i? s ), 
one can derive in case of the conical volume : 



w(0 = 



' con \ 
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2/3 
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(14) 



where the expression in backets is the argument with 
which W c ^i of Eq. flU) has to be evaluated. The lower in- 
tegral boundar y is now given by the value £ m in for which 
the relation £\/£ 2/3 - 1 = 1/3 is fulfilled. 

The relativistic effects listed above can be taken into 
account in the following approximate way. General rela- 
tivistic ray bending (point (3) in the list at the beginning 
of this section) requires to replace the Euclidian values 
of Ho and /ii of Eq. (|l^) by their relativistic counterparts. 
These are found by tracing geodesies from an observer po- 
sition on the axis to the points where they hit the torus 
surface tangentially. Figure |3^ shows two examples where 
we calculated the ray trajectories in the gravitational field 
of a black hole with approximately 3M Q for a torus with 
radius R t = 1R S and a torus with radius Rt = 3i? s , re- 
spectively. In our analysis we allow for only one torus 
image and thus neglect that neutrino trajectories may 
be bent around the black hole to cross the system axis 
more than once. Blueshift, redshift and radiation trap- 
ping (points (1), (2) and (4), respectively) are accounted 
for by assuming that they can be split off in a separate 
factor (S GR C GR ) in the integrand of Eq. (|ll|). Thus for 
the cylindrical integration volume we take 



w(0 = w^(0 = w^(0{s^-c^)(0, 



(15) 
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and for the integration over the cone we use 
W(0 EE w^(0 = 



(16) 

The factor (5 GR C GR ) in Eq. (|l5|) depends on £ and in 
Eq. ( |l6| ) on 3£-\/£ 2 ^ 3 — 1: and the brackets (...) indicate 
that the product S C has been averaged over the 
plane at constant z. For the latter averaging procedure 
one has to notice that S GR C GR , as specified below, has a 
weakly diverging singularity at the event horizon, but its 
integral has a finite value and is regular. 

The red- and blueshift factor S GR can be approxi- 
mated by 



5 GR (r) = e 2 ^ 



1 - 



e 0(r o ) 



Tf \5/2 / p x -3/2 
its \ ^ JXs 



3R S + R* 



1 - — 
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(17) 



where r denotes the radius where neutrinos and antineu- 
trinos annihilate, the lapse function is given by = 
-\/l — Rs/r, and the second expression is obtained by tak- 
ing r « 3i? s + Rt as the mean radius from where neutri- 
nos and antineutrinos originate. The first factor e 2 ^ * 1 in 
Eq. ([It]) accounts for the redshift and time dilation from 
r Q to infinity, and the second factor (e^ r °) /e^* 1 ) 3 for the 
blueshift or redshift between r Q and r. 

Close to the event horizon only photons moving nearly 
radially outward can escape from the gravitational attrac- 
tion of the black hole whereas photons propagating with 
large angles relative to the outward direction are captured 
by the hole. We assume that the electron-positron-photon 
plasma is locally isotropic and neglect hydrodynamic ef- 
fects. This allows us to estimate the escape probability 
of energy produced by vv annihilation at radius r from 
the fraction of photons on escape trajectories according 
to (Shapiro & Teukolsky |1983[) : 



C7 GR (r) 



1 




(18) 



with the minus sign holding for R s < r < | R s 
plus sign for r > |i? s . At r = |i? s a fraction of 50% of the 
pair plasma is able to escape to infinity. The assumption of 
local isotropy probably implies an underestimation of the 
escape probability because the non-zero radial momentum 
of pairs produced by vv annihilation is not taken into 
account. Ray bending towards the black hole and positive 
pressure gradients, on the other hand, act in the opposite 
direction and reduce the expansion of the pair plasma out 
of the close vicinity of the event horizon. Therefore we 
consider Eq. ( |l8| ) as a reasonable zeroth-order estimate. 

The expressions of Eqs. @ and (|l|) are now inserted 
into Eqs. (O) and (|l^) to evaluate the integral of Eq. (|TT|) . 
Equation (0) then yields the energy deposition rate by vv 



annihilation as measured by an observer at rest at infinity, 
when the neutrino luminosities and the mean neutrino en- 
ergies are taken as measured at the neutrino source. The 
annihilation luminosities including general relativistic cor- 
rections are reduced relative to the Newtonian results. The 
calculated reduction factors are listed in Table || for dif- 
ferent values of the parameter R t /R s . Relativistic effects 
decrease the available energy in all investigated cases by 
typically 10-30%, mainly because of redshift to infinity 
and radiation capture by the black hole. However, these 
two effects are nearly compensated by the increase of the 
vv annihilation probability due to the facts that the neu- 
trinos emitted towards the black hole are blueshifted and 
the vv annihilation rate rises with the square of the neu- 
trino luminosity times the average neutrino energy. It is 
interesting to note that general relativistic ray bending 
reduces the angle under which the torus surface is seen 
from a position on the system axis near the black hole, 
whereas the viewing angle of the torus is increased rel- 
ative to the Euclidean angle at large distances. For this 
reason, ray bending implies a reduction of the vv energy 
deposition rate in the case of the cylindrical integration 
volume, but has the opposite effect for the conical volume 
where a larger fraction of the total energy comes from vv 
annihilation far away from the black hole. 

7.3. Additional heating by neutrino- electron scattering 

The pair-plasma cloud which is created around the ac- 
cretion torus by vv annihilation may receive additional 
energy input by neutrino scattering on the abundant elec- 
trons and positrons (Woosley 1993a ). In this subsection we 
estimate the importance of this heating process relative to 
vv annihilation with and without relativistic corrections. 
We adopt again the simplified torus-black hole geometry 



described in Sect. 7.2 



In the pair-plasma cloud electrons and positrons are 
present in equal numbers and thus the electron chemical 
potential is zero. For such conditions the summed energy 
transfer rate by scattering of neutrinos and antineutrinos 
of all flavors becomes (Tubbs and Schramm 1975) 



and the Qe 
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with Vj — v e , v e , ZAj, v^, v T , v T . Here St is the surface area 
of the torus, er = 1-76 x 10~ 44 cm 2 , and c eUj are weak 
interaction coefficients. In order to estimate the electron 
energy density e e and the average electron energy (e e ) ~ 
4T, we assume that the e ± 7 plasma stays near the torus 
for an expansion time t cxp . Then the plasma temperature 
can be estimated by equating the energy density of the 
plasma with the heating rate times the expansion time. 
One gets 
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Using this in Eq. (^) and integrating over the volume 
outside of t he torus with the same assumptions as made 
in Sect. 7.2, we find 



^ £ 0.075 Lv " 62{e 



v. /10 



u 



1 ms ) a(3V ' 



(21) 



where the numerical value was obtained by taking rep- 
resentative numbers for the neutrino luminosities and 
mean neutrino energies from Figs. p3|-p6|. The normal- 
ized quantities are L v ^% — L„ c /10 52 ergs -1 , (e I/e )io = 
(e„ e )/10 MeV and R s g = i? s /9km. For the expansion 
timescale we assumed 1 ms which is probably an upper 
limit because the light crossing time through the main re- 
gion of neutrino energy deposition is only a few 10 s. 
The quantity V in the denominator is given by Eq. (pT|), 
and hi is defined by 
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d£(Mi-MoMO<^*)(0. 



(22) 



The factor (5 GR ) is a function of £ and different from 
unity only when general relativistic corrections are taken 
into account. It contains additional terms for the redshift 
between the radius of neutrino emission (r ), the radius 
of neutrino-electron scattering (r), and the observer at 
infinity. The brackets (...) indicate again that this factor 
is averaged on the planes of constant vertical height z. For 
the function 5 GR (r) we estimate 



5 GR (r) = e<« r °> 



1 - 



Re 



3R S + Rt 
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(23) 



In the general relativistic case, the time t exp is assumed 
to be measured by the observer at rest at infinity. 

In Table | numbers for the factor U/(a(3V) of Eq. (|l]) 
are listed for both the cylindrical and conical integra- 
tion volumes and with and without relativistic corrections. 
Since this factor is always much less than unity, the con- 
tribution to the pair-plasma heating by neutrino-electron 
and -positron scattering is unimportant relative to vv an- 
nihilation. Certainly it contributes not more than about 
10% of the e^7 energy for all parameter combinations sug- 
gested by our models. Again, general relativistic effects 
lead only to a moderate reduction of the results. 

8. Summary and discussion 

8.1. Summary 

With three-dimensional hydrodynamic simulations we 
continued our evolution calculations of merging neutron 

1997| , and in particular Ruf- 



stars (Ruffert et al. L996 



fert and Janka 1998b) into the phase when the massive 
(~ 3Mq baryonic mass) object at the center of the merger 



remnant has collapsed to a black hole, and the surround- 
ing material with high angular momentum is forming an 
accretion disk around this black hole. From the available 
set of merger models we chose a case where a massive 
disk could be expected, and varied the time after the 
coalescence when the black hole formation was assumed 
to happen. Moreover, we studied the subsequent accre- 
tion phase with two different prescriptions for the gravi- 
tational potential of the black hole, either of Newtonian 
or of Paczyhski-Wiita type. The latter potential is deeper 
and reproduces the existence of the innermost stable circu- 
lar orbit of a Schwarzschild black hole at 3R S . With these 
three-dimensional simulations we were not able to follow 
the long-time evolution of the accretion disk and the ac- 
cretion of its matter until completion, a process which is 
governed by the viscous transport of angular momentum. 
Nevertheless, we could study the approach to a nearly 
quasi-stationary state where the radial velocities are much 
smaller than the orbital velocities. 

We determined the disk properties in the quasi-sta- 
tionary state, such as its mass, temperature, density, an- 
gular momentum, neutrino emission, and rate of mass loss 
into the black hole and thus its lifetime. The numerical re- 
sults presented in this paper are in reasonably good agree- 
ment with analytical estimates on grounds of a very sim- 
ple torus model by Ruffert et al. (1997), and with recent 
calculations of hyper-accreting black holes by Popham et 
al. ( ^998| ). 

The disk properties and the computed disk evolution 
do not depend very much on the exact moment when the 
black hole formation is assumed to take place, unless this 
happens earlier than about 2 ms, or a very long time af- 
ter the two neutron stars have merged. We found, how- 
ever, that the amount of gas which does not fall into the 
black hole immediately on a dynamical timescale, but or- 
bits around the hole for at least several revolutions, de- 
pends strongly on the employed gravitational potential of 
the black hole. Part of this sensitivity is certainly due to 
the fact that our simulations of the neutron star coales- 
cence had been carried out with Newtonian gravity and 
therefore the use of the Paczyhski-Wiita potential for the 
torus models is not fully consistent. 

A reliable prediction of the mass which stays in the 
accretion torus for several orbital periods is possible by 
the criterion that the specific angular momentum of the 
matter must exceed the Keplerian value at radius 3R S . 
This means that the orbital velocities in the disk are 
only slightly sub-Keplerian. Pressure effects are not too 
important to determine the radial disk structure, but 
of course are crucial for the vertical disk structure. In 
case of the Paczyhski-Wiita potential a gas mass of only 
3.5 x 10 -2 M Q is able to complete more than 5 orbits 
around the black hole. A quasi-stationary state was not 
reached within the 5 ms of simulated disk evolution, and 
the high accretion rate of 7M Q s _1 implies an estimated 
accretion timescale of at least another 5 ms at the end 
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of the simulation. On the other hand, the final state of 
the disk with Newtonian potential after 5 ms shows very 
small radial velocities. The accretion rate is 5 M Q s -1 and 
the accretion of the disk mass of 2.6 x 1O -1 M0 can be 
estimated to continue for more than 50 ms (Table |||) . 

The quasi-stationary evolution of the torus is deter- 
mined by numerical viscosity. For the chosen grid res- 
olution and the input physics implemented in our hy- 
drodynamics code we estimate a viscosity parameter of 
a ~ 0.01. Maximum temperatures in the tori are around 
10 MeV, and maximum densities are between several 10 11 
gem -3 and 10 12 gcm -3 with somewhat higher values in 
the Newtonian case. Due to the larger mass, the Newto- 
nian torus is optically thick to neutrinos whereas in the 
Paczyhski-Wiita calculation the disk is nearly transparent 
for v e and heavy-lepton neutrinos. For this reason the to- 
tal neutrino luminosities in both cases are similar, around 
10 53 ergs -1 . The v e luminosity is about twice as high as 
the v e luminosity, and heavy-lepton neutrinos contribute 
only a minor fraction to the total energy loss. The mean 
energies of the emitted neutrinos are around 10 MeV for 
electron neutrinos, 13-16 MeV for electron antineutrinos 
and 15-20 MeV for the heavy-lepton neutrinos (Table |l|). 

Because of the rapid radial infall of the torus gas in 
the Paczyhski-Wiita simulation on the one hand, and the 
rather large neutrino opacity of the Newtonian torus on 
the other, the efficiency q u for transforming rest-mass en- 
ergy into neutrino emission is below the maximum theo- 
retical limit of 8.3% for the radiation efficiency of a New- 
tonian accretion disk (5.7% for relativistic disk accretion 
onto a nonrotating black hole). We found values for q v of 
0.5% for the Paczyhski-Wiita model and of 1.3% for the 
Newtonian computation (Table ||) . Neutrino-antineutrino 
annihilation in the surroundings of the black hole deposits 
energy at a rate of (3-5) x 10 50 ergs -1 in the region where 
the baryonic mass density is below 10 11 gem -3 . This cor- 
responds to an efficiency of 0.4-0.5% for the conversion 
of neutrino energy into e + e - pairs. A fraction of 10-30% 
of this vv annihilation energy or an estimated time inte- 
gral during the accretion phase of E v p ss 10 48 -10 49 erg are 
dumped in a cone around the system axis where the total 
baryon mass is below 10 -5 M Q . A comparison of Figs. [l7] 
and |3l] shows that the radius of this low-density cone at 
its base is 15-20 km. 

General relativistic effects like gravitational redshift 
and blueshift, radiation capture and ray bending reduce 
the observable energy in the expanding pair-plasma only 
insignificantly mainly because the different effects act in 
opposite directions and thus partly cancel each other. 
However, because the black hole has an initial relativis- 
tic rotation parameter a = Jc/(GMg H ) ~ 0.4 and is spun 
up to a w 0.5 by the accretion of the torus material, effects 
associated with the Kerr character of the black hole may 
lead to an increase of the lifetime of the torus and may 
raise the energy dumped in e + e - pairs above the values 
obtained in this work. This was discussed by Jaroszyhski 



( |1996| ) and Popham et al. ( 1998 ), although the results in 
these papers cannot be directly applied to the very high 
mass accretion rates of our models where the tori start to 
become opaque against neutrinos. 

The high baryon density and integral baryon mass out- 
side of an axis-near cone suggest that relativistic expan- 
sion of an e ± 7 plasma can only develop in jets along the 
axis. From the shape of the low-density funnel above the 
poles of the black hole (approximately bounded by the 
contour of density 10 9 g cm -3 in Figs. [l5| and |l^ or by the 
contour corresponding to 10 -5 Mq in the plots of the z- 
d-integrals of Fig. [l7]) we expect that the relativistic pair 
plasma will break out with rather wide opening half-angle 
between about 10 degrees and several 10 degrees, which 
implies a moderate beaming factor for the two jets of 
fa = 2SQ/{ATr) = 1 - cos0 between 1/10 and 1/100. This 
means that the determined vv annihilation energies may 
be sufficient to explain gamma-ray bursts with observed 
energies of E 1 w A-kE vV /{28Vl) < 10 50 -10 51 erg, if the ob- 
server assumes isotropy of the source. For short 7 bursts 
with typical durations i 7 between 0.1s and Is, luminosi- 
ties L 1 Ki E^/try between 10 50 ergs -1 and 10 52 ergs -1 may 
be within the reach of our models. If Kerr effects play a 
role even larger energies and luminosities may be possible. 

8.2. Discussion of implications and outlook 

Figure H|| sketches the envisioned geometry of the for- 
mation of the pair-plasma jets and their properties as 
suggested by our models. Energy deposition by electron 
neutrino and antineutrino absorptions in the outer lay- 
ers of the accretion disk is likely to drive a baryonic 
mass flow off the disk surface similar to the neutrino- 
driven winds of newly-formed hot neutron stars (Woosley 
& Baron |1992|, Qian & Woosley |1996|). This wind ex- 



pands with subrelativistic velocities and has much lower 
entropies than the 7-jets in which the entropy is ex- 
tremely high: sjet « (e + P)/(Tn b ) ~ f m u £U/(TM je t) « 
|(m u c 2 /T)r ~ 1000 T fc B per nucleon (e, P, n b , T are en- 
ergy density, pressure, baryon density and temperature, 
respectively, at the base of the jet, m u is the atomic mass 
unit, and Y the Lorentz factor of the jet). In and near 
the boundary layers between the wind region and the 
jets (hatched in Fig. |3^) a significant amount of dilute 
gas (~ 10 -3 Mq) may be present with such high temper- 
atures that the entropies could be much larger than in 
the neutrino-driven winds of nascent neutron stars. These 
ejecta from hot accretion tori around stellar mass black 
holes might therefore be an ideal site for high-entropy r- 
processing. Future simulations will have to determine the 
mass of this material and its contribution to the galactic 
nucleosynthesis of heavy elements. 

Despite of the promising results for energy and lumi- 
nosity of the gamma jets, our models still suffer from an 
unpleasant problem. The maximum relativistic Lorentz 
factors r — 1 = Evv/(Mc 2 ), which can be estimated from 
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the integral energy deposition rates given in Fig. ^ (mul- 
tiplied by the torus lifetime) and the integral masses of 
Fig. are only around 0.1-0.5. Even if the total energy 
deposited by vv annihilation gets somehow focussed into 
the jets, this estimate for T — 1 increases only to values 
between 1 and 5. This is roughly two orders of magnitude 
below the desired Lorentz factors of at least 100 which 
would require a baryon loading of at most 10 -7 Mq within 
the low-density cone around the system axis. The mass of 
roughly 10 -5 Mq that fills the region near the axis in our 
models, however, is only one order of magnitude above the 
numerical lower limit of the mass resolution of our simu- 
lations. The mass resolution is limited because we set a 
minimum baryonic mass density of p m i n = 5x 10 7 gcm -3 
in the cells of our grid. Therefore we cannot exclude that 
the large mass of dilute gas in the surroundings of our tori 
has purely numerical reasons. As a consequence, the polar 
regions of the black hole might not clear up in our models 
as fast as they actually would if the baryons were allowed 
to deplete unlimited because they are sucked in rapidly 
along the axis by the black hole. 

There is another, even more important effect which 
can help opening up a clean funnel for relativistic plasma 
jets. In a short but extremely luminous outburst of neu- 
trinos right after the neutron star merging, peak rates 
of the energy deposition by vv annihilation (and addi- 
tional v e and v e absorption) of more than 10 52 ergs -1 are 
reached. This produces very high energy densities above 
the polar caps of the compact massive object at the center 
of the merger remnant (for details, see Ruffert & Janka 
1998b). The maximum energy deposition rates are sev- 
eral 10 32 erg cm _3 s _1 at heights \z\ ps 30 km in gas with 
densities between 10 9 g cm -3 and 10 10 g cm -3 . This means 
that the equivalent of the gravitational binding energy is 
transferred to the matter within only one millisecond. The 
heated plasma must expand extremely rapidly, because 
it experiences only little resistance from the gas farther 
out due to the very steep density gradient and the corre- 
spondingly small density scale height above the poles of 
the forming black hole. The expanding hot gas will push 
overlying baryons away and the baryon den sities must 
drop quickly. Recently MacFadyen & Woosley ( 1998 ) have 
performed simulations of collapsing, rapidly rotating and 
massive stellar cores ( "collapsars" ) . Although the situa- 
tion there is less favorable because of the huge mass that 
surrounds the black hole and the accretion torus, their 
simulations demonstrate impressingly the depletion of the 
axis-near region due to the outward expansion of the gas 
in response to powerful energy input. This creates a much 
cleaner funnel for the pair-plasma jet which is driven by 
the annihilation of neutrinos and antineutrinos emitted 
from the disk around the black hole during the subse- 
quent, longer period of accretion. For all these reasons we 
do not consider the current problem of baryon pollution in 
the jets as seriously worrying as yet. This problem is most 
likely surmountable by simulations which achieve a better 



mass resolution and take into account the feedback effects 
of the vv annihilation around the torus self-consistently. 



Figure 34 provides an overview over different possible 



evolution paths of merging binary neutron stars and neu- 
tron star black hole binaries. Whether these systems lead 
to short gamma-ray bursts by vv annihilation or not de- 
pends on the parameters of the components (masses and 
spins) and on the properties of the nuclear equation of 
state. All events are accompanied by th e emission o f grav- 
itational w aves ( see, e.g., Shib ata et al 
& Shapiro |l994|; Zhuge et al 



1995 



1992 



1993 



Rasio 



Ruffert et al. g_996|) 



and neutrinos (Ruffert et al. 1997 ), and by the ejection of 
mass with possible implicat ions for heavy-element nucle- 
osynthesis (Rosswog et al. 1998 and references therein). 
Although Kerr effects associated with the rotation of the 
black hole can increase the energies relative to our results, 
vv annihilation is most likely not able to yield enough en- 
ergy for long complex gamma-ray bursts. Magnetic fields 
and magnetohydrodynamic conversion of gravitational en- 
ergy into radiation may be required to explain the long 
and complex gamma bursts by NS-NS and NS-BH merg- 
ers (Thompson 1994, Meszaros & Rees 1997, Meszaros et 
al. 1998), unless this subclass of bursts is connected with 



events described by the "failed supernova" or "collapsar" 
models (Woosley 1993a, Popham et al. 1998, Paczyhski 



1998, MacFadyen & Woosley 1998) 



An option also not included in the scheme of Fig. |34] is 
the one that the neutron stars collapse to black holes prior 
to their me rging . This scenario was sug gested by Mathews 
& Wilson (|l998j ) and Mathews et al. fll998j ) on grounds 
of results of general relativistic simulations which yield a 
compression of the stars as they spiral in towards each 
other, instead of tidal stretching as suggested by Newto- 
nian and post-Newtonian simulations. The possibility of 
a collapse depends on the question how close the neutron 
star is below the maximum stable mass for the consid- 
ered nuclear equation of state. It is currently a matter of 
vivid discussions whether the numerical results by Math- 
ews & Wilson ( 1998| ) and Mathews et al. (1998) are correct 
and in agreement with analytical understanding of the ef- 
fects of general relativity on the pre-merging evolution (see 
Thorne 



1997 and references therein). 



Similarly extreme phases of mass accretion by a black 
hole as studied in this paper occur during the merging of 
neutron star black hole binaries. Three-dimensional hy- 
drodynamic simulations of the dynamical interaction of 
these binary stars and of the post-merging evolution were 
performed by Eberl (1998a) with the same input physics 
as described in this paper and in Ruffert et al. ( |1996| , |1997D , 
i.e., gravitational- wave and neutrino emission were taken 
into account. The results will be published separately 
(Eberl et al. 1998b). They reveal that the accretion disks 
can be more massive (^ 0.5 M©) than those obtained in 
neutron star neutron star mergers. Since the average ac- 
cretion rates are significantly higher, the neutrino lumi- 
nosities are up to 10 times larger, but the accretion phase 
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is shorter, so that effectively the gamma jets receive about 
10 times more energy. The neutrino luminosities show con- 
siderable variation because of large temporal fluctuations 
of the accretion rate. 

Kluzniak & Lee fll998|) and Lee & Kluzniak ( |l998|) 
have discovered the possibility of perhaps long periods of 
episodic mass transfer from the orbiting neutron star to 
the black hole. These are connected with cycles of orbital 
decay due to gravitational-wave emission and subsequent 
widening of the orbit after mass has been donated to the 
black hole. The existence of these cycles might be sensitive 
to the properties of the nuclear equation of state (assumed 
to be a polytropic law in the simulations by Kluzniak & 
Lee 1998 and Lee & Kluzniak 1998 ) and on the binary 
parameters like the neutron star spin and the mass ratio 
of the neutron star to the black hole. 

If neutron stars near the minimum stable mass were 
formed after dozens or hundreds of these cycles (Klu zniak 
fc Le e 1998 ) and then exp lode (Colpi & Shapiro 



1991, 



1993 



198S 



Sumiyoshi et al. 1998 ), it is very unlikely that 
gamma-ray bursts are produced because large amounts of 
baryonic matter will be ejected into the surrounding space. 
Even if the mass donating neutron star were heated by vis- 
cous dissipation due to the action of tidal forces during the 
cyclic phases of decreasing and increasing orbital distance, 
there is hardly a chance to make gamma-ray bursts by the 
associated neutrino emission. The heating of the neutron 
star surface layers by neutrino absorption would drive a 
slow baryonic wind rather than create a relativistically 
expanding pair-plasma with a low enough baryon loading 
(Paczyhski 1990 and, in particular, Woosley 1993b). 



The same argument also applies if the stiffness of the 
supranuclear equation of state prevented the rapidly rotat- 
ing, hot remnant of the binary neutron star merger from 
collapsing to a black hole on the dynamical timescale. In 
that case the remnant would cool by emitting its huge 
gravitational binding energy in neutrinos, which are able 
to escape from the dense object only on the diffusion 
timescale of several seconds. Again, the surroundings of 
the merger would be polluted with the baryons which are 
driven off the surface of the remnant in a continuous flow. 

In contrast, the black hole-accretion disk systems pro- 
vide ideal conditions for efficient vv annihilation and the 



creation of relativistic o utflow (Woosley 1993a , Mochko- 
vitch et al. 1993, 1995). On the one hand, the neutrino 



source is very compact which ensures large neutrino num- 
ber densities. On the other hand, the axis-near region 
above the poles of the black hole is an environment with 
naturally low baryon loading where neutrinos and antineu- 
trinos can collide at large angles and annihilate with high 
efficiency. Our simulations have confirmed that accreting 
black holes from coalescing binary neutron stars can in- 
deed power short (~ 0.1s) gamma-ray bursts with lumi- 
nosities up to 10 52 ergs _1 , provided the jets are beamed 
into 1/100-1/10 of the sky. Our models need to be im- 
proved and extended with respect to general relativistic ef- 



fects, in particular associated with the possible Kerr char- 
acter of the black hole. In addition, it would be interesting 
to study cases with different disk viscosities and to include 
the effects of vv energy deposition self-consistently in the 
hydrodynamic calculations. 
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Fig. 15. Vertical structure of the Newtonian Model B10 in the x-z-plane and y-z-plane near the end of the simulation 
(t » 14.5 ms). The plots correspond to the cuts in the equatorial plane given in Fig. 11. The upper left figure shows 
contours of the baryonic mass density (measured in gem -3 ) together with the flow field indicated by velocity vectors, 
the upper right figure displays contours of the temperature (in MeV), the lower left figure of the electron fraction Y e , 
and the lower right figure of the entropy per nucleon. The density contours are spaced logarithmically with intervals 
of 0.5 dex (bold: log(p [gem -3 ]) = 9.0, 10.0, 11.0, 12.0), the temperature contours are spaced linearly, starting with 
1 MeV and 2 MeV, and then continuing in steps of 2 MeV (bold: T = 10 MeV), the Y e contours correspond to the 
levels 0.02, 0.04, 0.06, 0.08, 0.10, 0.15, 0.20, 0.25, and 0.30, and the entropy levels were chosen to be 1, 2, 3, 4, 5, 
6, 8, 10, 12, 14, 16, and 20 (bold: 6, 10, 20). Some of the contours are labeled with their respective values, and darker 
grey shading indicates higher values of a quantity. 
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Fig. 16. Vertical structure of the Paczyriski-Wiita Model B10 in the x-z-plane and y-z-plane near the end of the 
simulation (t 14.9 ms). The plots correspond to the cuts in the equatorial plane given in Fig. [l^. The upper left figure 
shows contours of the baryonic mass density (measured in gem -3 ) together with the flow field indicated by velocity 
vectors, the upper right figure displays contours of the temperature (in MeV), the lower left figure of the electron 
fraction Y e , and the lower right figure of the entropy per nucleon. The density contours are spaced logarithmically with 
intervals of 0.25 dex (bold: log(p [gem -3 ]) = 9.0, 10.0, 11.0), the temperature contours are spaced linearly, starting 
with 1 MeV and 2 MeV, and then continuing in steps of 2 MeV (bold: T = 10 MeV), the Y e contours correspond to the 
levels 0.02, 0.04, 0.06, 0.08, 0.10, 0.15, 0.20, 0.25, and 0.30, and the entropy levels were chosen to be 1, 2, 3, 4, 5, 
6, 8, 10, 12, 14, 16, and 20 (bold: 6, 10, 20). Some of the contours are labeled with their respective values, and darker 
grey shading indicates higher values of a quantity. 
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Fig. 17. Mass distribution in the Model B64 of Ruffert & Janka ( 1998b ) at t = 10 ms (upper figure), in the Newtonian 
Model B10 at time t = 14.9 ms (lower left figure) and in the Paczyhski-Wiita Model B10 at time t = 15.2 ms (lower right 
figure). In the upper panels of the figures ("z-integral"), the azimuthally averaged baryon densities were integrated 
along the z-direction from a given value of z out to the grid boundary (practically infinity because of the rapidly 
decreasing density) and within radial intervals of width Ad — 1 km. In the lower panels ("z-d-integral"), integration in 
radial direction from to the value of d on the abscissa was additionally performed. The baryonic mass is measured 
in Mq and the corresponding contours of equal-mass levels are spaced logarithmically in steps of — 1. Therefore each 
point (rf, z) in the upper panels is associated with the mass inside a hollow cylinder between d — Ad and d from z 
to infinity, and in the lower panels with the mass inside a full cylinder between d = and d from z to infinity. The 
contours are mirrored along the d — axis in all plots and show only the mass in the volume on one side of the 
equatorial plane. The evacuated funnel along the system axis is clearly visible in the final stages of Models B10 and 
B10. 
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Fig. 18. The dots represent the azimuthal velocities v v (d) 
(normalized to the speed of light) of all zones in the equa- 
torial plane of the neutron star merger Model B64 of 
Janka & Ruffert (1998b) at t = 9.34ms after the start 
of that simulation. This model was close to the one (at 
t = 10.0 ms) used as initial condition for the tori evo- 
lution simulations, Models B10 and B10. The thin solid 
line gives the average value of the azimuthal velocities 
(binned in steps of 3 km), and the bold solid line the lo- 
cal Newtonian Kepler velocity fK ep ier(^) as a function of 
the equatorial radius d. 

azimutal velocities, model B10 
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Fig. 20. The dots represent the azimuthal velocities v v (d) 
(normalized to the speed of light) of all zones in the equa- 
torial plane of the Paczyhski-Wiita torus Model B10 at 
t = 14.94 ms near the end of the simulation. The thin 
solid line gives the average value of the azimuthal ve- 
locities and the bold solid line the local Kepler velocity 
^KcpicrW °f the Paczyhski-Wiita potential as a function 
of the equatorial radius d. The kinks of the curves inside 
the inner vacuum boundary at two Schwarzschild radii 
(2R S ) of the central black hole (about 18 km) are unphys- 
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Fig. 19. Same as Fig. 18 but for the Newtonian torus 
Model B10 at t = 14.51 ms near the end of the simulation. 
The kinks of the curves inside the inner vacuum boundary 
at two Schwarzschild radii of the central black hole (about 
18 km) are unphysical and caused by the interpolation of 
the plot data. 
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Fig. 21. Average specific angular momentum j(d) — 
v ip (d)d (in units of km times the speed of light) of the gas 
in the equatorial plane of the Newtonian torus Model B10 
(bold solid line) at t = 14.51 ms and of the Paczyhski- 
Wiita torus Model #10 at t = 14.94 ms (thin solid line) 
as function of the equatorial radius d. The specific angular 
momentum jK ep i er (d) corresponding to the Kepler veloc- 
ity of the Newtonian gravitational potential is displayed 
as bold dotted line, the equivalent quantity JKepier^) 
for the Paczyhski-Wiita potential as thin dotted line. 
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Fig. 22. The dots give the radial velocities v r (d) (normalized to the speed of light) as a function of the equatorial radius 
d for all grid zones in the equatorial plane of the Newtonian Model B10 at t = 14.51 ms (left) and of the Paczyhski- 
Wiita torus Model B10 at time t = 14.94ms (right). The solid lines represent the mean values of all zones within 
binning intervals of 3 km. In the Newtonian model the velocities are very small for d <; 25 km with a radial average of 
about — 10 cms , indicating that the model has relaxed to a quasi-stationary state near the end of the simulation. 
In the Paczyhski-Wiita torus the gas collapses rapidly inside the last stable circular orbit at d = 3R S « 27 km, but 
expands in the region 40 km <J d <J 80 km due to ongoing viscous heating and outward transport of angular momentum. 
Note that for d J> 100 km the gas is very dilute whereas the Newtonian torus extends beyond d = 160 km (see Figs. [l5| 
and 16). In both models the plotted values at d < 20 km are not physical due to the interpolation of the data. 
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Fig. 23. Luminosities of electron neutrinos (labeled with Fig. 24. Same as Fig. £3] but for the Paczyhski-Wiita 
e), electron antineutrinos (labeled with a), sum of all Models B2 and BIO. 
heavy- lepton neutrinos (labeled with x), and total neu- 
trino luminosity (labeled with tot) as functions of time 
for the Newtonian Models B4 and BIO. 
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Fig. 25. Average energies of the emitted electron neu- 
trinos (labeled with nu), electron antineutrinos (anti nu), 
and heavy-lepton neutrinos (mu, tau) as functions of time 
for the Newtonian Models B4 and B10. 



Fig. 26. Same as Fig. g5| but for the Paczyhski-Wiita 
Models B2 and BIO. 
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Fig. 27. Surface emissivity of energy in neutrinos (measured in ergcm _2 s _1 ) for Model BIO at the end of the simulation 
(time in the top right corner of the panels). The plotted contours were obtained by integration of the local energy 
loss rates from z = to infinity, in panel a for electron neutrinos, in panel b for electron antincutrinos, in panel c for 
the sum of all heavy- lepton neutrinos, and in panel d for the summed contributions from all flavors of neutrinos and 
antincutrinos. The contours are logarithmically spaced with intervals of 0.5 dex, bold contours are labeled with their 
respective values. The grey shading emphasizes the emission levels, dark grey corresponding to the strongest energy 
loss by neutrino emission. The dotted line marks the intersection of the z = plane with the surface of optical depth 
unity where neutrinos start to stream off freely. 
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Fig. 28. Same as Fig. gfl but for Model BIO at time t = 15.26 ms. 
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Fig. 29. Contour plots of the local energy loss rates (in ergcm~ 3 s _1 due to the emission of neutrinos and antineutrinos 
of all flavors. The left plot shows cuts in the x-z- and y-z-planes perpendicular to the equatorial plane of the Newtonian 
Model B10 near the end of the simulation, the right plot the corresponding cuts for the Paczyhski-Wiita Model ,810. 
The contours are spaced logarithmically with intervals of 0.5 dex, bold contours correspond to the values 28, 30 and 
32. The grey shading emphasizes the emission levels, dark grey corresponding to the strongest energy loss by neutrino 
emission. 
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Fig. 30. Maps of the local energy deposition rates (in erg cm -3 s -1 ) by vv annihilation into e + e~ pairs in the sur- 
roundings of the accretion torus of the Newtonian Model B10 at time t — 14.94 ms (left) and of the Paczyhski-Wiita 
Model B10 at time t = 15.26ms (right). The values shown as solid contour lines in a plane perpendicular to the 
equatorial plane were obtained as averages over the azimuthal angle around the z-axis. The abscissa d measures the 
distance from the grid center in the cc-y-plane. The white octagonal area around the center with a semidiamcter of one 
Schwarzschild radius indicates the presence of the central black hole. The contours are logarithmically spaced in steps 
of 0.5 dex, the grey shading emphasizes the levels with dark grey meaning high energy deposition rate. The energy 
deposition rate was evaluated only in that region around the torus where the baryon mass density is below 10 11 g cm -3 . 
The integral value of the energy deposition rate at the displayed time is 4.9 x 10 50 ergs _1 for Model B10 (left) and 
3.1 x 10 50 ergs _1 for Model B10 (right). The dashed lines mark the (approximate) positions of the neutrino "spheres" 
of v e , v ei and v x (from outside inward), defined by the requirement that the optical depths in z-direction are T ZVi = 1 
(see also Eq. (Q) and Figs. 27 and ^). The dotted contour lines represent levels of constant values of the azimuthally 
averaged mass density, also logarithmically spaced with intervals of 0.5 dex, the bold dotted line corresponding to 
p = 10 10 gcm" 3 . 
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Fig. 32. Geodesies of (massless) neutrinos emitted tangentially from a toroidal accretion disk (white circle) with 
assumed semidiameter i?t of one Schwarzschild radius, R s — 9 km, (left) and of three Schwarzschild radii (right) in 
the vicinity of a black hole (black circle) with a mass of 3M©. The torus center is at x = 3R S + Rt, the black circle 
has a radius of IR S . 
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Fig. 33. Envisioned gamma-jets and neutrino-driven baryonic wind from the massive accretion torus around the BH. 
The dashed circle indicates a region around the black hole with the radius of the last stable circular orbit. The hatched 
areas represent schematically the volume where baryons might be able to penetrate into the jet region; its outer 
boundary is approximately described by Eq. (|l0|), its inner boundary may roughly be given by the line where the 
gravitational force is balanced by the component of the centrifugal force in the opposite direction, fa = 2<50/(47r) 
denotes the fraction of the sky into which both jets are beamed. 
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Fig. 34. Possible evolution paths from NS-NS and NS-BH mergers to short gamma-ray bursts powered by neutrino- 
antincutrino annihilation into electron-positron pairs. The evolution of the merging binaries depends on the parameters 
of its stellar components (masses, spins) and on the properties of the nuclear equation of state. 



